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Abstract: 

We review aspects of loop quantum gravity in a pedagogical manner, with the aim of 
enabling a precise but critical assessment of its achievements so far. We emphasise 
that the off-shell ('strong') closure of the constraint algebra is a crucial test of quantum 
space-time covariance, and thereby of the consistency, of the theory. Special attention is 
paid to the appearance of a large number of ambiguities, in particular in the formulation 
of the Hamiltonian constraint. Developing suitable approximation methods to establish 
a connection with classical gravity on the one hand, and with the physics of elementary 
particles on the other, remains a major challenge. 
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1 Key questions 

When four-dimensional Einstein gravity is quantised canonically, in a perturbation series in the 
Newton constant around flat space-time, divergences arise at two-loop order. An impressive cal- 
culation by Goroff and Sagnotti P, 01 has demonstrated that in order to obtain a finite S-matrix, 
the action should contain a counterterm 

^'^^ = 72880 (16^2)2 I^^^VaC^^paC'^^^Cxr'"' , (1.1) 

a result which was later confirmed in an independent background-field calculation by van de 
Ven 01 • The usual conclusion drawn from this result, and from the fact that the coupling constant 
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is dimensionful, is that an infinite number of counterterms is needed. It is generally agreed that 
this non-renormalisability renders perturbatively quantised Einstein gravity meaningless as a 
fundamental theory because an infinite number of parameters would be required to make any 
physical prediction. 

However, this still leaves open the possibility that Einstein gravity can be quantised consis- 
tently, but that it is simply the perturbation series in Newton's constant which is ill-defined. 
This possibility has been raised and advocated not only in the context of canonical quantisation, 
but also, and independently, in the context of suggestions that there may exist a non-trivial fixed 
point of the renormalisation group in Einstein's theorv ^ [Hj 0] • Implicitly, it also underlies the 
path integral approach to Euclidean quantum gravity [3, la] ) which provides a possible framework 
for the discussion of semi-classical states jSl]. Indeed, to date there is no proof that such a quan- 
tisation which does not make use of a series expansion around a fixed background is guaranteed 
to fail. However, quantising without relying on perturbation theory around a free theory is hard. 

Loop quantum gravity, or LQG for short, is an attempt to quantise Einstein gravity non- 
perturbatively. In contrast to string theory, which posits that the Einstein-Hilbert action is 
only an effective low energy approximation to some other, more fundamental, underlying theory, 
LQG takes Einstein's theory in four spacetime dimensions as the basic starting point. It is a 
Hamiltonian approach, which is background independent in the sense that its basic quantities and 
concepts do not presuppose the existence of a given spatial background metric. In comparison to 
the older geometrodynamics approach (which is also formally background independent) it makes 
use of many new conceptual and technical ingredients. A key role is played by the Ashtekar 
variables, which allow the reformulation of gravity in terms of connections and holonomies, and 
which (at least in the form originally proposed) greatly simplify the constraints. A related key 
feature is the use of spin networks, which in turn requires other mathematical ingredients, such as 
non-separable ('polymer') Hilbert spaces and representations of operators, which are not weakly 
continuous. Undoubtedly, novel concepts and ingredients such as these will be necessary in order 
to circumvent the problems of perturbatively quantised gravity, but, as for any other approach 
to quantum gravity, it is important not to loose track of the physical questions that one is trying 
to answer. 

The goal of the present paper is to review some essential properties of loop quantum gravity 
in an easily accessible way from a non-specialist's perspective, and with a non-LQG audience in 
mind. It is not meant to be a comprehensive review of the subject; readers who want to know 
more about the very latest developments in this field are instead referred to a number of excellent 
recent reviews describing the story from the specialist's point of view jiol. 
we would like to provide an entree for 'outsiders', and to focus on the outstanding problems as 
we perceive them, and thereby initiate and enable an informed debate between the proponents of 
the different approaches. Accordingly, we will take the liberty to omit mathematical details that 
in our opinion are not truly essential to understand the physical consequences of the formalism 
(at least not on a first perusal of the literature), and to describe some results 'our own way'. 
At the same time, as we move along, we will try to make precise and clearly state the questions 
that are often raised about the LQG programme (for earlier reviews which mention some of these 
concerns, see 

In order to focus the discussion, and for the reader's convenience, we begin with a summary 
of what we consider to be the main issues and open questions. 

• How do Einstein's equations appear in the classical limit? 

The space of quantum states used in LQG (not necessarily a Hilbert space) is very different 
from the one used in Fock space quantisation. As a consequence, it becomes non-trivial 
to see how semi-classical 'coherent' states can be constructed, and how a smooth classical 
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spacetime might emerge. In simple toy examples, such as the harmonic oscillator, it has been 
shown that the LQG quantisation method indeed leads to quantum states whose properties 
are close to those of the usual Fock space coherent states In full (3+l)-dimensional 
LQG, the classical limit is, however, far from understood (so far only kinematical coherent 
states are known 0, El, El, 13 ) • In particular, we do not know how to describe or 
approximate classical spacetimes in this framework that 'look' like, say, Minkowski space, 
or how to properly derive the classical Einstein equations and their quantum corrections. 
A proper understanding of the semi-classical limit is also indispensable to clarify the con- 
nection (or lack thereof) between conventional perturbation theory in terms of Feynman 
diagrams, and the non-perturbative quantisation proposed by LQG. 

• Renormalisation vs. regularisation: where is the 2-loop divergence? 

No approach to quantum gravity can claim success that does not explain the ultimate fate of 
the two-loop divergence In a consistent scheme, this divergence must be either elim- 

inated by cancellation, or disposed of by some other mechanism. The key question raised 
by (|1.H1 is therefore what happens when one expands the results of non-perturbatively 
quantised Einstein gravity in Newton's constant. When such an expansion is performed 
about a semiclassical state (which remains to be found in LQG, see above), the two-loop 
divergence should manifest itself in one form or another. In its present incarnation, LQG 
cannot (yet?) 'see' and accommodate this divergence. Its possible 'disappearance' is occa- 
sionally argued to be due to the emergence of an effective cut-off (regulator) which might 
eventually make the perturbation theory finite. If this is the case, an obvious question 
concerns the true origin of this cut-off: is it generated dynamically, or covertly put in 'by 
hand' (for instance, by working with the compact group SU(2) instead of the full Lorentz 
group)? There are also questions concerning the meaning of 'regularisation'. According to 
conventional (quantum field theory) wisdom, physics is supposed not to depend on the way 
in which the theory is regulated before the cutoff is removed; how can it be that physics 
predictions of LQG do depend on the chosen regularisation prescription? This question is 
in part answered by the fact that the notions of 'finiteness' and 'regulator independence' 
as currently used in LQG on the one hand, and in conventional quantum field theory and 
perturbative quantum gravity on the other, are not the same; see section 1^31 Let us fur- 
thermore stress that in spite of its perturbative origin, the result cannot be so easily 
dismissed as a background artifact: while it does require some background for its derivation 
(i.e. a spacetime solving Einstein's equations), the counterterm is actually independent of 
the particular background about which one expands, see 0], as is also evident from the 
manifestly space-time covariant form in which it is written. 

• Status of the Hamiltonian constraint? 

In the current LQG literature, there is surprisingly little discussion of certain basic aspects 
concerning the Hamiltonian constraint operator, which are of central importance for the 
theory (recent exceptions are [2^, and the so-called 'master constraint programme' of 13). 
For this reason, we will here describe the Hamiltonian constraint operator and its action on 
a given spin network wave function in rather pedestrian detail, as far as we have been able to 
work it out. In particular, we will exhibit the numerous choices and ambiguities inherent in 
this construction, as well as the extraordinary complexity of the resulting expression for the 
constraint operator. The number of ambiguities can be reduced by invoking independence 
of the spatial background and indeed, without making such choices, one would not even 
obtain sensible expressions, as we shall see very explicitly. In other words, the formalism 
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is partly 'on-shell' in that the very existence of the (unregulated) Hamiltonian constraint 
operator depends very delicately on its 'diffeomorphism covariance', and the choice of a 
proper 'habitat', on which it is supposed to act in a well defined manner. A further source 
of ambiguities, which, for all we know, has not been considered in the literature so far, 
consists in possible /i-dependent 'higher order' modifications of the Hamiltonian (|5.14j) . 
which might still be compatible with all consistency requirements of LQG. 

The attitude often expressed with regard to the remaining ambiguities is that they corre- 
spond to different physics, and therefore the choice of the correct Hamiltonian is ultimately 
a matter of physics (experiment?), and not mathematics. We disagree, because we cannot 
believe that Nature will allow such a great degree of arbitrariness at its most fundamental 
level: recall that it was precisely the infinite number of parameters and the concomitant 
ambiguities which killed perturbative quantum gravity! 



Does the quantum theory possess full spacetime covariance? 

Spacetime covariance is a central property of Einstein's theory. Although the Hamiltonian 
formulation is not manifestly covariant, full covariance is still present in the classical theory, 
albeit in a hidden form, via the classical (Poisson or Dirac) algebra of constraints acting 
on phase space. However, this is not necessarily so for the quantised theory. LQG treats 
the diffeomorphism constraint and the Hamiltonian constraint in a very different manner. 
Why and how then should one expect such a theory to recover full spacetime (as opposed to 
purely spatial) covariance? The crucial issue here is clearly what LQG has to say about the 
quantum algebra of constraints. Unfortunately, to the best of our knowledge, the 'off-shell' 
calculation of the commutator of two Hamiltonian constraints in LQG - with an explicit 
operatorial expression as the final result - has never been fully carried out. Instead, a 
survey of the possible terms arising in this computation has led to the conclusion that the 



commutator vanishes on a certain restricted 'habitat' of states j27l. l28l. l29l| . and therefore 
the LQG constraint algebra closes without anomalies. By contrast, we will here argue that 
this 'on shell closure' is not sufficient for a full proof of quantum spacetime covariance, but 
that a proper theory of quantum gravity requires a constraint algebra that closes 'off shell', 
i.e. without prior imposition of a subset of the constraints. The fallacies that may ensue 
if one does not insist on off-shell closure can be illustrated with simple examples. In our 
opinion, this requirement may well provide the acid test on which any proposed theory of 
canonical quantum gravity will stand or fail. 



• Matter couplings: anything goes? 

Because LQG is claimed to be a finite and fully consistent theory of quantum gravity, it does 
not appear to impose any restrictions on the types of matter that are coupled to gravity, 
nor on their interactions. Indeed it is straightforward, though sometimes cumbersome, to 
extend the formalism to include matter: in this perspective, matter appears to be a mere 
accessory that can be added on to pure gravity as one chooses. This is in marked contrast to 
supergravity and superstring theory, which are based on the hypothesis that the very raison 
d'etre of matter is its indispensability for curing the perturbative (and non-perturbative) 
inconsistencies of quantum gravity, and the desire to 'geometrise' matter in the framework 
of a (probably supersymmetric) unified theory. It is not inconceivable that LQG might 
eventually encounter new consistency requirements when descending from the kinematical 
to the physical Hilbert space, but we see no evidence for this so far. Therefore the question 
remains whether and how LQG can recover the consistency requirements that conventional 
perturbative quantum field theory imposes on the matter content of the world, in particular 
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those resulting from cancellation of local (gauge) anomalies. Let us recall that Nature does 
'care' about such consistency requirements, in that it has chosen to put the three known 
generations of fermions into anomaly free multiplets of the standard model gauge group. 
Similar comments apply to global anomalies. Is it possible to obtain the correct answer for 
pion decay when fermions are coupled to electromagnetism in the LQG approach, or would 
LQG predict the neutral vr meson to be a stable particle? 



Structure of space(-time) at the smallest scales? 

There is a general expectation (not only in the LQG community) that at the very shortest 
distances, the smooth geometry of Einstein's theory will be replaced by some quantum 
space or spacetime, and hence the continuum will be replaced by some 'discretuum'. LQG 
does not do away with conventional spacetime concepts, in that it still relies on a spatial 
continuum as its 'substrate'. At the kinematical level, it imposes a discrete structure which 
is very different from the discreteness of a lattice or naive discretisation of space (i.e. of a 
finite or countable set), by 'polymerising' the continuum via the scalar product 1)4. 7|) . This 
is similar to the discrete topology ('pulverisation') of the real line with countable unions 
of points as the open sets. Because the only notion of 'closeness' between two points in 
this topology is whether or not they are coincident, whence any function is continuous in 
this topology, this raises the question as to how one can recover conventional notions of 
continuity in this scheme. 

However, the truly relevant question here concerns the structure (and definition!) of physical 
space and time. This, and not the kinematical 'discretuum' on which holonomies and spin 
networks 'float', is the arena where one should try to recover familiar and well-established 
concepts like the Wilsonian renormalisation group, with its continuous 'flows'. Because the 
measurement of lengths and distances ultimately requires an operational definition in terms 
of appropriate matter fields and states obeying the physical state constraints, 'dynamical' 
discreteness is expected to manifest itself in the spectra of the relevant physical observables 
(the associated kinematical operators will be discussed in section HT!]) . It is not clear whether 
discreteness of (physical) space would entail a discrete structure for time also, hence space- 
time, as there is no a priori notion of 'time' in quantum gravity. Instead, 'time' must 
also be defined operationally in terms of a 'clock field', see e.g. jiol . 31, [s^]. Continuity 
or discreteness of physical space and time thus follows from the properties of the relevant 
'measuring rod fields' and 'clock fields', and their spectra. 



• Conceptual issues 

Last but not least: although we will have nothing new to say here on the grand conceptual 
issues of quantum gravity and quantum cosmology we wish to emphasise that these issues 
must be addressed and resolved by all approaches to quantum gravity. This comment 
concerns not only the difficult interpretational problems but also more technical aspects. 
Among the former, we would like to mention the problem of interpreting the 'wave function 
of the universe' and associated 'matrix elements' between different such wave functions, 
or the problem of constructing 'observables' and their physical interpretation; among the 
latter, there is the question of whether we have any right to expect the 'wave function of the 
universe' to be normalisable, or the related question of whether the familiar Hilbert space 
formalism of standard quantum mechanics is really the correct mathematical framework 
for quantum gravity.^ Over the past decades, there has been slow but steady progress on 



In ordinary quantum mechanics, the Hilbert space formahsm and the concomitant notion of 'unitarity' are 
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several fronts (see e.g. 0, [s^, 31, [s^, 0, EB, E^, E^l and references therein), but it is 
probably fair to say that we are still far from fully understanding these issues. 



2 Topics omitted 

As we mentioned already, there are several recent developments and advances which we cannot 
cover here for lack of space, or lack of expertise on our side. We will list these and briefl y co mment 
on them below, but otherwise must refer readers to the pertinent 'inside' reviews jiol . 111 . l3 . 1^ . 
[l3 | and the more recent original references for more information and other points of view. 



2.1 Loop Quantum Cosmology 

Symmetry reduced versions of LQG have recently been studied as models of quantum cosmol- 
ogy ji^ . Is^. Two main features stand out. The first is a new mechanism for the avoidance of 
the big bang singularity in the framework of mini-superspace models of quantum gravity.^ This 
mechanism is based on the fact that the inverse scale factor is represented by an operator that 
stays bounded as the classical radius of the universe shrinks to zero. Alternatively, one might say 
that the effective discretisation of the Hamiltonian constraint enables the quantum wave function 
to 'jump over' the singularity, and to continue to evolve past the singularity. However, it is not 
clear whether and how these models can be derived from full-fledged LQG, and whether singu- 
larity avoidance is also a property of the full theory. In fact, a very recent investigation 41] has 
revealed that the spectrum of the operator corresponding to the inverse volume is not bounded 
from above in full LQG. In addition, it is by no means excluded (and some would even say likely) 
that inhomogeneities, and possibly other degrees of freedom, will essentially alter the nature of 
the quantum state near the singularity. See also for other comments. 

The second new feature arises in applications of loop quantum cosmology to inflationary mod- 
els. It is the possibility that inflation might be triggered and eventually stopped (gracefully) by 
gravity itself, via an intrinsically quantum gravitational mechanism |43|. Even if scalar fields 
cannot be avoided altogether, an appealing feature here is that the inflaton potential engineer- 
ing characteristic of most current models of inflation could become unnecessary. Perhaps less 
attractive (to an outsider) is that the regularisation ambiguities of LQG feed through to the sym- 
metry reduced models and lead to physical effects (see |4J| for a recent discussion). Concerning 
the value of the Barbero-Immirzi parameter, obtained from black hole entropy calculations (see 
below), one might arrive at an interesting 'internal' consistency check by matching it with the 
duration of the inflationary period, which is also linked to this parameter in these models. It has 
furthermore been suggested that measurements of the CMB fluctuation spectrum at large angles 
might provide experimental tests of LQG, but it is not clear to us whether the predicted effects, 
if present, might not be explained in many other ways, too. 



2.2 Microscopic origin of black hole entropy 



The explanation of the Bekenstein-Hawking entropy of black holes l4f)l. I4f 
has been claimed as a success for both LQG 0,0 

and string theory [49| 

of string theory is that it not only explains the area law, but also predicts the factor 1/4 relating 



in terms of microstates 
The main achievement 



tightly linked to the physical interpretation of the theory in terms of probabilities and their conservation in time. 
In the absence of an a priori notion of time (as in quantum gravity) it is therefore by no means evident whether 
these will remain the relevant concepts. 

^Let us note that 'singularity avoidance mechanisms' may exist also in more conventional mini-superspace 
quantum geometrodynamics . 
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entropy and area; furthermore, the agreement has been shown to extend to the higher order 
curvature terms predicted by string theory, see e.g. [HO]. However, the argument requires a huge 
extrapolation in the string coupUng constant, and is essentiahy Umited to BPS-type extremal 
black holes. The LQG explanation, on the other hand, succeeds in implementing the condition 
for isolated horizons [HH in the quantum theory, and works for ordinary (Schwarzschild and Kerr) 
black holes, but requires a 'fit' for the Barbero-Immirzi parameter to get the prefactor right. 
Although the two ansatze thus both reproduce the desired result, we are faced with something 
of a paradox here, as the two explanations seem almost impossible to reconcile, given the very 
different hypotheses underlying them — pure gravity on the one hand, and an exponentially 
growing spectrum of D-brane states on the other hand. 

2.3 Spin foams 

Attempts to overcome the difficulties with the Hamiltonian constraint have led to yet another 
development, spin foam models [HI, IH^, EJ]. These can be regarded as space-time versions of 
spin networks, to wit, evolutions of spin networks in 'time'. Mathematically, these models rep- 
resent a generalisation of spin networks, in the sense that group theoretical objects (holonomies, 
representations, intertwiners, etc.) are attached not only to vertices and edges (links), but also 
to higher dimensional faces in a simplicial decomposition of space-time. In addition, they may 
pave the way towards the construction of a physical inner product, which in turn is defined by 
the quantum dynamics. Many of the recent advances in this area concern purely topological 
theories, so-called models", where F{A) is a field strength, and B the Lagrange multiplier 
(tensor) field whose variation enforces F(A) = 0. The formalism can thus be nicely applied to 
(2+1) gravity, which is topological. In (3+1) dimensions one needs extra constraints in order to 
recover the propagating gravitational degrees of freedom back into the theory [H^. Interesting 
as they are, however, these developments have so far not shed much new light on the problems 
with the Hamiltonian constraint, or the constraint algebra, in our opinion. A derivation of spin 
foam models from the Hamiltonian formulation remains incomplete due to the complexity of the 
Hamiltonian constraint [s^. Hence, a decisive proof of the connection between spin foam models 
and the full Einstein theory and its canonical formulation still appears to be lacking. On the 
contrary, it has even been suggested that these models may provide a possible 'way out' if the 
difficulties with the 'orthodox' Hamiltonian approach, which we follow here, should really prove 
insurmountable. 



3 Old vs. new variables: from metric to loops 
3.1 Prelude: the metric (or dreibein) approach 

There are many introductory texts to which we can refer readers for a more detailed treatment of 
canonical gravity (see e.g. 3, |^, El, El, E^, E^, ) , but let us nevertheless briefly review the 



traditional way of doing canonical quantum gravity, also called geometrodynamics. Our exposition 
in this section follows whose notations and conventions we adopt in the remainder. We will 
be using the vierbein formalism with a vierbein -E"^^ and space-time metric G^i, = E^'^Ey^rjAB 
with flat (tangent space) metric tiab = diag (— 1, +1, +1, +1) ^. The vierbein is covariantly 
constant under a derivative which is covariant w.r.t. both spacetime diffeomorphisms and local 
Lorentz transformations, viz. 

D^E,^ := V^^/ - LO^^^'E^B = , (3.1) 



Modulo some dimension dependent factors, the results described in this subsection are valid in any dimension. 
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where is the covariant derivative which involves only the Christoffel connection, 



V^K := - r^/l^p (3.2) 

and ujf_iAB is the spin connection. It is a standard result that both T^j^'' and uj^ab can be explicitly 
solved for the vielbein from the above equation (in the absence of torsion). 

As is customary in canonical gravity, and following the standard ADM prescription [H^l we 
assume the space-time manifold M to be globally hyperbolic. It then follows from a Theorem of 
Geroch [^^l that can be foliated according to = S x M, with a spatial manifold S of fixed 
topology (and no boundary, for simplicity). Using letters m, n, ... and a, 6, ... for curved and flat 
spatial indices, respectively, we choose the triangular gauge for the vierbein 

V^(^ fl) ^ G,.-(--^-^--" ;»), (3.3) 

where Qmn = ^m^^na is the spatial metric. The vierbein components N and Na are Lagrange 
multipliers ('lapse' and 'shift'); their variation yields the canonical constraints. In the remainder 
we will freely convert spatial indices by means of the spatial dreibein and its inverse. 

The canonical momenta are obtained in the standard fashion from the Einstein-Hilbert action 

= ^eer (i^af, - W , (3.4) 

where e = det(em"), and 

Kab '■= Ea^Efy^V fj^E^o = iOabO (3.5) 

is the extrinsic curvature of S in 7W (expressed in terms of flat spatial indices, with ujabo = 
Ea^cOnbo), and K = Kaa- The inverse formula reads 

Kab = -{Tiab-\5ab^) , (H = H,,) . (3.6) 
From the symmetry Kab = we immediately deduce the Lorentz constraint 

Lab = em[a^b]"^ ~ , (3.7) 

which is also the canonical generator of spatial rotations on the dreibein (« means 'weakly 
zero' [^). The canonical Hamiltonian (really Hamiltonian density) is 

H = dtea'^Um" -C = NHo + NaHa , (3.8) 

with the diffeomorphism constraint 

Ha = DraUa"^ - , (3.9) 

and the Hamiltonian constraint {alias the scalar constraint) 

Ho = {Uabllab - ^tf) - eii(3) = ^e{KabKab - K^) - ei2(3) ~ q , (3.10) 

where i?*-'^^ is the spatial Ricci scalar. The canonical equal time (Poisson) brackets are 

{e^a(x),nfe"(y)} = ,5,6C,5(3)(x,y) , (3.11) 
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with the other brackets vanishing (for the canonical variables with the indices in the indicated 
positions). Canonical quantisation in the 'position space representation' now proceeds by repre- 
senting the dreibein as a multiplication operator, and the canonical momentum by the functional 
differential operator 

na'"(x) = --4^. (3.12) 

With these replacements, the classical constraints are converted to quantum constraint operators 
which act on suitable wave functionals. The diffeomorphism and Lorentz constraints become 

//„(x)^[e]=0 , L„f,(x)^'[e] =0. (3.13) 

They will be referred to as 'kinematical constraints' throughout. Dynamics is generated via the 
Hamiltonian constraint, the Wheeler- DeWitt (WDW) equation 6^, 6^, 6^, 6^ : 



i?o(x)*[e] =0. (3.14) 

It is straightforward to include matter degrees of freedom, in which case the constraint operators 
and the wave functional \I'[e, . . . ] depend on further variables (indicated by dots). The functional 
^'[e, . . . ] is sometimes referred to as the 'wave function of the universe', and is supposed to contain 
the complete information about the universe 'from beginning to end'. A good way to visualise ^ 
is to think of it as a film reel; 'time' and the illusion that 'something happens' emerge only when 
the film is played. 

The substitution (|3.12|) turns the Hamiltonian constraint into a highly singular functional 
differential equation, which most likely cannot be made mathematically well defined in this form, 
even allowing for certain 'renormalisations'. Here we do not wish to belabour the well known 
difficulties, both mathematical and conceptual, which have stymied progress with the WDW 
equation for over forty years, see e.g. [3, |3) ^- We should point out the (surely well 
known) fact that, at least formally and ignoring subtle points of functional analysis, solutions 
to both the Lorentz constraint and the diffeomorphism constraint can be rather easily obtained 
by integrating suitably densitised invariant combinations of the spatial dreibein (or metric) and 
curvature components, and the matter 'position variables' over the spatial three-manifold S. 
The associated wave functionals ^ are then automatically functionals of diffeomorphism classes 
of spatial metrics. In saying this, it must be stressed, however, that geometrodynamics has so 
far not succeeded in constructing a suitable scalar product and an appropriate Hilbert spaces of 
wave functionals. 

The absence of a suitable Hilbert space is often invoked by LQG proponents as an argu- 
ment against the geometrodynamics approach. While LQG is certainly 'ahead' in that it does 
succeed in constructing a Hilbert space (more on this later), we should like to emphasise that 
in all approaches there looms the larger conceptual problem of whether conventional quantum 
mechanical concepts are indeed sufficient for quantum gravity (31i. i33i]. For instance, even if we 
can compute matrix elements of wave function(al)s, we still have no idea what their correct phys- 
ical interpretation is in the context of quantum cosmology, or whether the normalisability of the 
wave function of the universe is a necessary condition. In other words, we do not know whether 
these cherished concepts may not have to be amended or abandoned in the final theory. For this 
reason, we believe that besides the emphasis on mathematical rigour, it is equally important to 
develop some physical intuition for the states one is dealing with, and in this regard, we do not 
think that geometrodynamics lags behind. From this point of view, it appears to us that, beyond 
the technical subtleties, the kinematical constraints are not the real problem of quantum gravity. 
The core difficulties of canonical quantum gravity are all connected in one way or another to the 
Hamiltonian constraint - irrespective of which canonical variables are used. 



*One of the inventors of 13.141 has been overheard aptly referring to it as "that damned equation.. 
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3.2 Ashtekar's new variables 

Much of the initial excitement over Ashtekar's discovery [b^ of new canonical variables was due 
to the change of perspective they bring about, which fuelled hopes that they might alleviate some 
of the longstanding unsolved problems of quantum gravity. Let us therefore first describe what 
they are, and how they are obtained. There is an ab initio derivation based on the addition of a 
term ocjEAEARto the Einstein-Hilbert action Ig^ . |^ (this term vanishes upon use of the 
Bianchi identity and the equations of motion for the spin connection with vanishing torsion), but 
we will skip this step here. Instead consider the combination ^ 

-^ma '■— ~ 2^abc^mbc ~l" 7^ O^ma ~ '^^ma^) — ~ '^^abc^mbc "l~ T-^ma ; (3.15) 

where uJmbc is the spatial spin connection. T he p arameter 7 7^ is referred to as the 'Barbero- 
Immirzi parameter' in the LQG literature T^, 3- Classically, 7 has no physical significance, but 
is believed to become physically relevant upon quantisation, for instance, bv setting the scale for 
the fundamental areas and volumes (in this sense it is somewhat analogous [T^] to the 6 parameter 
of QCD). One can now show that j67i | 

{Ama{^), Anbiy)} =0, 

{^r(x), K(y)} = 0, (3.16) 

{A„'^(x),K(y)} = 7W('^(x,y), 

where the canonically conjugate variable to A^a is the inverse densitised spatial dreibein 

E^:=eea'^ =^ E = e\ (3.17) 

with e = detem" and E = detE'^. The parameter 7 is often eliminated from these brackets by 
absorbing it into the definition of E^. 

To rewrite the constraints in terms of the new variables, we first observe that the covariant 
constancy of the spatial dreibein and the Lorentz constraint imply the Gauss constraint: 

DmK = dmK + eabcAjE'^'^ « . (3.18) 

Defining the field strength 

Pmna ■ — "^m^na f^n^ma ~l~ ^abcAmbAnc 

1 \ 2 b c (^-^^^ 

— '^^abcRmnbc ~l" ^(-^m-^na I^nKrna) ~\~ 7 ^abc-^m I^n j 

it follows that the diffeomorphism constraint takes the form 

E:^F^na ~ . (3.20) 

Furthermore, 

eabcETKFmnc = -^^UabHab " ^tf) - R^^^ 



2 



(3.21) 



^We alert readers that our notational conventions differ from the ones used in most of the LQG literature, 
where a,b, . . . denote curved and i, j, . . . flat spatial indices. 
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These relations immediately suggest interpreting as a gauge connection for the gauge group 
SO (3) of spatial rotations (for D = 2 + 1 gravity, this group is replaced by its non-compact form 



S0(l,2) [73])- Accordingly, the new variables are conveniently rewritten as |67| 

= Ama'''af3 i (3.22) 

where are the Pauli matrices. 

For the special choice 7 = iti [^3] the second term on the r.h.s. of (|3.21|) drops out, and — 
save for an extra density factor of e — the Hamiltonian constraint is expressed entirely in terms 
of the new canonical variables, and furthermore depends on them polynomially.^ In other words, 
this particular choice allows us to combine the diffeomorphism and Hamiltonian constraints, 
which are schematically of the form '(5 + uj)Ii ~ 0' and '(5u; + lj^ + H^) ~ 0', respectively, 
into a single expression in terms of the connection A = uj + ill, and its canonically conjugate 
variable. Moreover, for this choice of 7 the connection Ami is nothing but the pullback of the 
four-dimensional spin connection to the spatial hypersurface S, with the two signs corresponding 
to the two chiralities (indicated by superscripts ^^^) 



'mABl^'' = u;„a67"' + 2^„^ao7'^7° = ^^1^7%"(1 + 7^) + ^L-i7°7'^(l - 7') , (3.23) 



(cf. eqn. (|3.5|) ). Equivalently, the variables Ami are associated with what, in a Euclidean formu- 
lation, would be the selfdual and anti-selfdual parts of the spin connection, respectively, 

^i^i = '^So with J^Ij, := ^ (^UJmAB ± '^eABCDU^m^''^ ■ (3.24) 

This is therefore also the natural choice for 7 when one considers coupling gravity to chiral 
fermions. In fact, one of the authors (H.N.) was first enticed into learning about Ashtekar's 
variables when realising that they simplify the calculation of the constraint algebra of supergravity 
considerably [2^,EHl- This is because the local supersymmetry constraint (i.e. the time component 
of the Rarita-Schwinger equation) always contains a factor D[miA)'^n]j where tpn is the gravitino, 
and A is just the Ashtekar connection (of course with 7 = iti). More succinctly, in supergravity, 
the commutation relations 1)3. 16() . as well as the polynomiality of the constraints must necessarily 
hold, if the commutator of two local supersymmetry constraints is to close into the scalar and 
diffeomorphism (and possibly other) constraints. 

From the esthetical point of view 7 = iti is therefore clearly the preferred choice. Nevertheless, 
this value has been abandoned in most of the recent LQG literature, because there is a major 
difficulty with it: the phase space of general relativity must be complexified with imaginary or 
complex 7. To recover the real phase space of general relativity and to ensure that real initial data 
evolve into real solutions, suitable reality conditions must be imposed. This is straightforward to 
achieve for the classical theory — after all, we have merely changed the variables, not the theory 
itself — but not so for the quantum theory. There, the complexification poses subtle problems 
concerning the definition and imposition of appropriate hermiticity conditions on the states and 
operators, and no consensus has been reached so far on how to circumvent these difficulties (or 
on whether they can be circumvented at all). 

For this reason one now usually takes 7 to be real. In this case, no problem arises with 
reality of solutions in either the classical or the quantum theory, but the polynomiality of the 
constraints, and hence one of the most attractive features of the new variables, is lost. This is 



^More recently it has, however, been appreciated that the extra density factor in 13.2111 is a major source of 
problems for a background-independent quantisation, since it makes the Hamiltonian an object of density weight 



two; for details see ,74J. 
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because the extra term in (|c{.21|) no longer vanishes, but must be subtracted from both sides to 
recover the correct WDW Hamiltonian. Accordingly, one must express the extra contributions in 
terms of the new canonical variables A^"' and E"^ via (|3.15|) and (|3.17|) . and this in turn requires 
expressing the original dreibein as well as the extrinsic curvatures in terms of the new canonical 
variables. For a while this was regarded as a chief obstacle, until a way to solve it was discovered 



by Thiemann 7J]. To this aim let us first introduce the volume associated with a region C S 



(considered as a phase space variable) ^ 



V{n) = jj?xe = jj?x \/I = ^d^x ^y^npe^bcErnEnEP . (3.25) 
Writing V = VCS), we first use the substitution 

e^'^(x) = emnpe^^'E-'/^Km^) = ^{^™"(x), v} (3.26) 

to recover the spatial dreibein. The second trick is to eliminate the extrinsic curvature using a 
doubly nested bracket. The first bracket is introduced by rewriting 

i^^'^(x) = -{A^"(x) , K} where K = K{^) := / d^xK^^Ea"" . (3.27) 
The second bracket comes in through identity 

= -17^1 e'^'^Fmnci^) , V} . (3.28) 



Last, we need another dreibein factor to convert the curved index m on to a flat one, for 
which we must use ()3.26() once again. We will postpone writing out the Hamiltonian constraint 
in terms of these multiple brackets until section \^7Il But let us note already here that the above 
relations supply the 'tool kit' also for transcribing any given matter Hamiltonian in terms of the 
new variables. 



For 7 7^ ibi the connection is no longer the pullback of a four-dimensional spin connection 77 1 . 
Although this does not immediately lead to problems with spacetime covariance, we will see in 
section 14.61 that the problem comes back when one considers coupling the theory to fermionic 
matter degrees of freedom. 

3.3 The connection representation 

Independently of the choice of 7, the reformulation of canonical gravity in terms of connection 
variables opens many new avenues, in particular the use of concepts, tools and techniques from 
Yang-Mills theory. Early attempts at quantisation were based on the connection representation 
with the original Ashtekar connection, i.e. j = ±i. Although these were ultimately not successful, 
let us nonetheless briefly summarise them here. In this scheme, one represents the connection 
Am"" by a multiplication operator, and sets 

The WDW functional depending on the spatial metric (or dreibein) is accordingly replaced by 
a functional ^[A] living on the space of connections (modulo gauge transformations). The price 



''The e-symbol is always the invariant tensor density: emnp = e ^em''eri''ep'^eabc, i.e. assumes the values 0,±1. 
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one pays is that this representation is much harder to 'visuahse' because the spatial metric is no 
longer represented by a simple multiplication operator, but must now be determined from the 
operator for the inverse densitised metric 



Even if one ignores the clash of functional differential operators at coincident points, finding 
suitable states and computing their expectation values is obviously not an easy task (and has not 
been accomplished so far). Similarly, the spatial volume density is obtained from 

,(.).^(x).!f.-.„,^^^ , ,3.31, 

Again this operator is very singular. Equations (|3..S0j) and (|3.3H) provide a first glimpse of the 
difficulties that the formalism has in finding semi-classical states and thereby establishing a link 
between the quantum theory and classical smooth spacetime geometry. 

For the quantum constraints the replacement of the metric by connection variables leads to 
a Hamiltonian which is simpler than the original WDW Hamiltonian, but still very singular. 
Allowing for an extra factor of e (and assuming e / 0, oo) the WDW equation becomes 



e 



Here we have adopted a particular ordering, which however is by no means singled out. No viable 
solutions to this constraint have been found, but there is at least one interesting solution if one 
allows for a non-vanishing cosmological constant A. Namely, using an ordering opposite to the 
one above, and including a term h.g with the volume density (|3.31|) . the WDW equation reads 

This is solved by 

^a[^] = exp J^d^^CcsiA)^ , (3.34) 

with the Chern-Simons Lagrangian £cs = AAdA+iAAAAA (actually, is already annihilated 
by the operator in parentheses, so the first factor in the Hamiltonian constraint operator (|3.33 | 1 
is not 'needed' for this result). In the literature this state is known as the Kodama state |78ii] . 
but the solution had been known for a long time in Yang-Mills theory [t^, where however it 
has rather unusual physical properties js^]- The difficulties with this solution have been much 
discussed recently |8ol . IsH : an obvious one concerns the flat space limit A ^ (idem for the 
'semiclassical' limit ^ ^ 0). 

What happens when we choose 7 to be different from zizi, and real in particular? As we 
explained already, the extra term in 1)3. 21() then no longer vanishes, must be dealt with sep- 
arately The nice polynomial form of the Hamiltonian constraint operator (|3.32j) is lost. 
When implementing the translation rules at the end of the foregoing subsection in the connection 
formulation, one finds that the new Hamiltonian is not significantly simpler any more than the 
original one of geometrodynamics in terms of metric or dreibein variables. 



14 



3.4 From connections to holonomies 

The loop representation is an attempt to overcome the difficulties with the connection representa- 
tion which we sketched above. The transition between the connection and the loop representation 
was originally obtained via the loop transform, which can be thought of as a kind of functional 
Fourier transform [s^l- We will not describe that construction here, but turn immediately to the 
formulation in terms of holonomies, on which the modern formulation of LQG - spin networks 
and spin foams - are based. 

Whereas in the connection representation one works with functionals ^[A] which are sup- 
ported 'on all of S', one now switches to the holonomies as the basic variables. These are gauge 
covariant functionals supported on one-dimensional links, or 'edges', which we will designate by e 
(following established LQG notation). For a given edge, i.e. some (open) curve embedded in S, 
we set 

he[A] =Vexp I Amdx'^ , with A = A^Ta . (3.35) 

Hence, he[-A\ is a matrix valued functional. The holonomy transforms under the action of SU(2) 
at each end of the edge e: 

hM] ^ K[A] = 5(e(0)) K[A] g-\e{l)) , with g{e{^)) , g{e{l)) G SU(2) . (3.36) 

For the remainder it is important that the holonomies are to be regarded as variables in their 
own right, subject to these transformation properties (so in some sense one can 'forget' about the 
original connection A defining the holonomy). The distributional nature of the holonomy is not 
only evident from its singular support (on a line rather than all of S), but also from the fact that 
we do not assume he.[A\ to be close to the identity if the edge e is 'small' (this terminology has to 
be used with due care, as there is no a priori measure that tells us when e is 'small', but we can 
still imagine making it 'small' by chopping e into as many 'subedges' as we like). The fact that 
the typical field configuration is generically a distribution rather than some smooth function is 
well known from constructive quantum field theory |^83]. 

The holonomies are taken to transform in SU(2) representations pj^ of arbitrary spin je = 
^, 1, |, . . . for each link e (with the convention that = means that there is no edge). We will 
denote such a spin-je valued holonomy by 

[p^Xhe[A\))^f,, (3.37) 

with indices a,(3, . . . as appropriate for the representation at hand. To make the notation less 
cumbersome, we will occasionally suppress pj^ and the representation labels, and simply denote 
the above matrix as (/ie[74])a^. 

To define the conjugate variable, we recall that the area element for the spatial manifold S 
can be expressed as a Lorentz vector {i.e. with flat S0(3) indices) via 

dFa := CabcO^ A 0^ with 0" = dx'^em" ■ (3.38) 

Happily, this can be nicely rewritten in terms of the new canonical variables 

dFa = e,fe,e„^e/dx" A dx^ = e^^pS^dx" A dx^ . (3.39) 

As the conjugate variable to h(,[A] one takes the 'flux' vector 

F§[E] := / dF^ (3.40) 
Js 
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Figure 1: Setup used for the computatfon of the bracket H3.42() . In the hmit e ^ the edge 
e(e) shrinks to zero and the two nodes just above and below the surface coincide. 



through any two-dimensional surface S embedded in E. There is also a smeared version of this 
variable, with a test function fa to soak up the free index a, which reads 

Fs[EJ]:= [ reabcen''e/dx^AdxP= [ /.e^^^^^Mx" A dx^ (3.41) 
Js Js 

We note that the standard notation for this variable in the LQG literature is E[S,f], but we 
prefer the one above because it is in parallel with the notation for the holonomy itself. Both 
Fs[E, f] and Fg[E] are distributional in the sense that they are supported on a two-dimensional 
submanifold of S. 

To compute the Poisson brackets between the new canonical variables introduced above, we 
consider a surface S and an edge e that 'pierces' S at the point P (if e does not intersect S, the 
bracket simply vanishes). We next subdivide this edge into three pieces as shown in figure ^ 
two subedges ei and 62, with associated holonomies and hl^, which touch S only in the limit 
e — > 0, and a third 'infinitesimal' edge e(e) intersecting S, for which the path ordered exponential 
can be approximated by the linear term. Then 



{ihM]U,Fs[E,f]} 



lim 



{KM])^A [ drr-A^<^(x)r^", , / dy" A dy%„p,/,(y)^''^(y) 1 (/i^JA]) 

I J e(e) J S I 



Sf3 



(3.42) 



L{e,Shfa{P){K,[A]T-heM])^p- 



Here, the intersection number 

i{e, S):= jdx^ dyM/ e„„p 5^^^ (x, y) = ±1 or , (3.43) 

encodes the information on how e intersects S* in a coordinate independent way. The integral 
is equal to ±1, depending on the orientation of e and 5", if e intersects S transversally. When 
the edge e intersects or touches S tangentially, a little care must be exercised; one then finds 
that t(e, 5) = 0, and the above bracket vanishes The fact that he[A\ and Fs[A^ f\ are, 

respectively, supported on one-dimensional and two-dimensional subsets of S is thus precisely 
what is required to perform the integral over the three-dimensional 5-function. Evidently, the 
integral is ill-defined when an entire segment of e lies within S] for a discussion of how to deal 
with this difficulty, see 
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4 Quantisation: kinematics 



Having determined the classical canonical variables, one would now like to promote them to quan- 
tum operators obeying the appropriate commutation relations. The essential assumption of LQG 
is that this quantisation should take place at the level of the bounded hermitean operators he[A\ 
rather than the connection A itself. This is analogous to ordinary quantum mechanics, when 
one replaces the Heisenberg operators x and p by Weyl operators e*^' and e*^; the spin network 
representation act ually uses the analog of a hybrid formulation with x and e*^. The Stone- von 
Neumann theorem 8^, 8^, S^l is usually invoked to argue that it makes no difference whether one 
quantises the Heisenberg or the Weyl algebra, i.e. that these quantisations are equivalent. The 
theorem does require, however, that the representations which are used are 'weakly continuous'. 
In the case of ordinary quantum mechanics, for example, this means that matrix elements of the 
operators corresponding to e*"^' and e^^^ are smooth functions of the parameters a and /?. In 
LQG the representations of operators do not satisfy this requirement.^ 

The failure of operators to be weakly continuous can, as we will see, be traced back to the very 
special choice of the scalar product ()4.7() below, which LQG employs to define its kinematical 
Hilbert space Wkin- This Hilbert space does not admit a countable basis, hence is non-separable, 
because the set of all spin network graphs in T, is uncountable, and non-coincident spin networks 
are orthogonal w.r.t. (|4.7|) . Therefore, any operation (such as a diffeomorphism) which moves 
around graphs continuously corresponds to an uncountable sequence of mutually orthogonal 
states in Wkin- That is, no matter how 'small' the deformation of the graph in S, the associated 
elements of Wkin always remain a finite distance apart, and consequently, the continuous motion 
in 'real space' gets mapped to a highly discontinuous one in "Hkin- Although unusual, and perhaps 
counter-intuitive, as they are, these properties constitute a cornerstone for the hopes that LQG 
can overcome the seemingly unsurmountable problems of conventional geometrodynamics: if the 
representations used in LQG were equivalent to the ones of geometrodynamics, there would be 
no reason to expect LQG not to end up in the same quandary. 

It is perhaps also instructive to contrast the LQG approach with the standard lattice approach 
to field theory. In the lattice approach, all quantities depend explicitly on the lattice spacing 
(i.e. the regulating parameter). In the limit in which the lattice spacing is taken to zero, one 
recovers the continuum results (all expectation values are smooth functions of the regulating 
parameter). In the LQG approach the 'discretuum' is instead built in by the very construction 
of the scalar product, rather than by introducing a regulating parameter. LQG is therefore quite 
different from 'quantum gravity on the lattice'. While the radical modification underlying LQG 
has certain appealing properties, it makes it hard to recover long-distance physics and the usual 
notion of continuity. 

In short, one thus implements the quantisation not by the replacement (|!'{.29|) but rather by 
promoting the above Poisson bracket ()3.42|) to a quantum commutator: 

\he[A]),,f,,Fs[E,f]\ =imhLie,S)riP){heAA]TaheM]),,p (4.1) 

or equivalently, 

\hM])ai3 , F§[E]\ = ihlh i{e, S) [K, [A]r" K, [A]) . (4.2) 

On the spin network representation the holonomies h(.[A\ will be represented as multiplication 
operators; the action of the canonically conjugate operators will be explained below. The in- 



*This is also the reason why the kinematical Hilbert space employed in loop quantum cosmology is different 
from the standard one already for a finite number of degrees of freedom When the number of degrees of 

freedom is infinite (as in quantum field theory), the Stone- von Neumann theorem anyhow does not apply. 
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equivalence of LQG quantisation with Fock-space quantisation arises through a special scalar 
product, to be discussed below. 

4.1 The Hilbert space of spin networks 

After defining the basic variables in which the theory should be quantised, the next step is to 
choose a Hilbert space in which the operators act. Starting from this space, one should construct 
a Hilbert space of physical states, i.e. space of states for which all constraints hold. The initial 
Hilbert space of LQG is the space of spin networks. While the Gauss constraint is easily solved in 
this space it turns out that a solution of the diffeomorphism constraints lies outside this 'naive' 
initial space, and one is forced to introduce a larger space. 

The intuitive idea behind spin networks is that the geometry at the Planck scale is foam-like. 
Physical gravitational and matter degrees degrees of freedom are excited only on so-called spin 
networks, i.e. one-dimensional 'edges' or 'links', and the vertices connecting them. Geometries 
which look smooth at large scales are supposed to arise only from complicated spin network 
states with many edges. In order to find the Hilbert space of these objects, one has to find a 
basis of wave functions over the configuration space, which associate a complex number to each 
and every configuration of the gauge connection. LQG makes use of wave functions which have 
singular support in the sense that they only probe the gauge connection on one-dimensional 
networks embedded in the three-dimensional spatial hypersurface S, which is a (not necessarily 
differentiable) manifold, i.e. can be mapped out by local charts. This three dimensional 'reference 
space', or 'carrier space' of the spin networks, does not carry any physical metric. LQG makes 
occasional use of local coordinates, or fiducial background metrics for certain intermediate steps 
in the construction, but physical quantities must not depend on such background data. We will 
here avoid their use as far as possible. Let us emphasise again that the 'discreteness' of the 
spin networks does not correspond to a naive discretisation of space. Rather, the underlying 
continuum, on which the spin networks 'float', the spatial manifold S, is still present. As we will 
see, the setup furthermore requires the a priori exclusion of infinite spin networks, that might 
contain Cantor-like or 'fractal' sets. 

By definition, each network is a (not necessarily connected) graph F embedded in S and 
consisting of finitely many edges Cj G F and vertices v G T. The edges are connected at the 
vertices. Each edge e carries a holonomy he[A] of the gauge connection A (this connection does 
not have to be smooth on the edge). The wave function on the spin network over the graph F 
can be written as 

^rAA]=^{heAA], K,[A],...), (4.3) 

where the tp is some function of the basic holonomies associated to the edges e G F. If, in addition, 
the wave function ■0 is invariant under arbitrary SU(2) gauge transformations it satisfies the Gauss 
constraint, and vice versa. A gauge invariant function ip thus takes care of joining the collection 
of holonomies into an SU(2)-invariant complex number by contracting all SU(2)-indices of the 
holonomies with invariant tensors 'located' at the vertices v, see figure |2 The basic building 
blocks of the spin network wave functions are therefore expressions of the following type. A 
three-valent vertex connects three edges according to 

ij[figm = (pniheAA])) (pniKM])) (pjAheA^])) , C^f.Sa • • • ' (^.4) 

V / ai/3i \ / a2P2 ^ ^ asps ^'-^^^•^ 

where dots represent the remainder of the graph. While the contraction is obviously unique if 
only two or three edges meet at a vertex, there may be more and independent choices for vertices 
of valence four and higher, depending on the way in which the edges are connected with the 
Clebsch-Gordan coefficients. For this reason, any given spin network will in general admit several 
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Figure 2: Examples of spin network states. For the 3-valent vertices on the left, the three 
incoming edges at each vertex are connected by a Clebsch-Gordan coefficient. For the 4-valent 
vertex on the right, one has to decide on a given way to construct a higher-order invariant tensor 
from two Clebsch-Gordan coefficients. 



independent wave functions of the above type. As an example, consider a four-valent vertex: one 
first has to decide on how to pair the edges into two groups of two. One such choice leads to e.g. 



\ / aifii \ / 02P2 ^ ^ 03P3 ^ ^ 04P4 

In these equations C'^^'^j^ Eire intertwiners (Clebsch-Gordan coefficients). In 1)4. 5|) the intermediate 
spin k can be freely chosen in accordance with the standard rules for the vector addition of angular 
momenta: |ji — J2I ^ k < ji +j2- In other words, we can graphically represent the 4-valent vertex 
by splitting it into two 'virtual' 3-valent vertices and adding a 'virtual' edge, carrying angular 
momentum k (see figure 13 • The same wave function can be re-expressed by performing this split 
in a different 'channel' by means of recoupling relations for the Clebsch Gordan coefficients (see 
for instance ;8Sj |) 

ci'''^ci'^f^ = yv2kTTV2^^l^} ^]ci'^'"ici'^r^, (4.6) 

0102P 0304P / ^ ^ 1 73 74 I "i'^2P 03a4P ' ^ '' 

where the object with curly brackets is the Wigner 6j-symbol. To make the dependence of a 
general spin network wave functional on the spins je associated with the edges and the intertwin- 
ers Cy associated with the vertices explicit, one occasionally designates functionals such as ()4.3|) 
also by |r, {j}, {C}). The space of finite linear combinations of such states is denoted by S (for 
'spin networks'). Although three- valent spin networks are obviously simplest, we will see later 
(see section 14.4.2(1 that higher valence is even generic because three- valent networks correspond 
to 'zero volume', and hence are deemed to be not of much interest. 

The wave functionals (|4.3|) are called cylindrical, because they probe the connection A only 
'on a set of measure zero' (like the (5- function does for ordinary functions). Similar 'cylindrical 
functionals' were formerly used in constructive quantum field theory in order to rigorously define 
the functional measures for free and certain interacting models of quantum field theory as limits 
of finite dimensional integration measures [s^] • The space S spanned by finite linear combinations 
of such cylindrical functions over all possible graphs is the starting point for the construction of 
the Hilbert space of spin networks. Obviously, the product of two cylindrical functions supported 
on the same or different spin networks is again a cylindrical function. 
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(a) 



(b) 



(c) 



Figure 3: A four-valent vertex is defined by a particular way of connecting tliree-valent vertices. 
Tlie spin k has to satisfy the triangle inequalities, but is otherwise arbitrary. 

To complete the definition of the space of spin network states, we must introduce a suitable 
scalar product. In LQG this is not the standard scalar product induced by a Fock space represen- 
tation (see the following section for more on this); instead, the scalar product of two cylindrical 
functions ^'r,{j},{c}[^] and ^r',{j'},{C'} [^] is defined as 

{^r,{j},{c} I ^r',{i'},{c'}) 

' if T^T', 

[ Yl<iK ^r,{j},{c}[he„...^i^r,^{fy^{C'}{hei,---) if T = F , ^^'''^ 

/ e,er 

where the integrals J dhg are to be performed with the SU(2) Haar measure. The above scalar 
product may look artificial, but given a few reasonable assumptions, there in fact exists a strong 
uniqueness theorem; the reader is referred to 0, E3, 0, S 0, for details. 

From ()4.7|1 we see that the inner product vanishes if the graphs F and F' do not coincide 
(even if they are 'very close to each other' in any given fiducial background metric). If F and 
F' coincide, the product may still vanish, depending on the choice of spins and intertwiners. At 
any rate, for F = F', the inner product is given by integrating the product of '0 and '0' over the 
holonomies at all edges, using the standard Haar measure. Because F was assumed to consist 
of finitely many edges, this is always a finite- dimensional integral, with one SU(2) integral for 
every edge in the graph. It is noteworthy that the above scalar product is invariant under spatial 
dijfeomorphisms, even if the states ^'i and ^'2 themselves are not, because the statement whether 
two graphs coincide or not is diffeomorphism invariant. Whereas the original wave functionals 
^[A] probe the value of the connection and therefore also depend on the position of the graph, this 
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information is 'lost' in the scalar product (|4.7jl which makes no more reference to the underlying 
space of connections or the 'shape' of the spin network graph. 

One now defines the kinematical Hilbert space Wkin as the completion of the space of gauge 
invariant spin network states S w.r.t. to this norm. Wkin thus consists of all linear superpositions 
of spin network states = |r„, {jn}, {Cn}) such that they have finite norm, 

oo 
n=l 

where the norm is defined by (|4.7|) . One of the distinctive features of the Hilbert space 

T^kin; which will play a key role in the further development of the theory, is its non-separability. 
This non-separability can be traced back to the existence of an underlying continuum, the spatial 
manifold S, and accounts for the difference between this scheme and the standard (Fock space) 
quantisation of gauge theories, see the following subsection. Although each in the above sum 
is associated to a finite graph r„ the expected number of edges need not be finite, because the 
sum 

oo 

(#(edges)) oc Y,K\^ L{Tnmn\W (4.9) 

n=l 

with L{Tn) the number of edges in r„, can be made to diverge if L{Tn) increases sufficiently 
rapidly with n, even if X]„ jonP I |^n| P < oo. Idem for the expectation value for the number 
of vertices. Let us, however, caution readers already at this point that Wkin is not the relevant 
Hilbert space for solving the quantum constraints in LQG. The physical Hilbert space, consisting 
of those states which satisfy the constraints, is expected to be separable; we will have more to 
say about this in the following sections. 

A second unusual feature, at least for a traditional quantum field theorist, is the ab initio 
absence of negative norm states, despite the fact that at this stage the constraints have not even 
appeared yet. This is in stark contrast to the usual covariant Fock space quantisation of gauge 
theories, where negative norm states are unavoidable, and can only be eliminated (if they can 
be eliminated at all) by restricting the Hilbert space to the subspace of physical states, which 
by definition are annihilated by the constraints. Curiously, there does exist a reformulation 
of (free) second quantised Maxwell theory with similar properties jo^. Roughly speaking, this 
consists in (i) defining the scalar product to be identical with the standard Fock space one for 
gauge invariant expectation values (as is well known, in second quantised Maxwell theory this 
restriction indeed eliminates negative norm states), and (ii) by declaring gauge- variant scalar 
products to be zero. In this way one preserves Lorentz covariance (because a gauge variant state 
is never transformed into a gauge invariant one by a Lorentz transformation, and vice versa), 
and eliminates negative norm states, in such a way that for observable quantities the results are 
identical with the usual Fock space quantisation. The price one pays is the loss of weak continuity 
(because the expectation values for gauge variant expressions always vanish, no matter how 
close they are to gauge invariant ones) and the loss of separability of the Hilbert space. These 
are precisely the features encountered in the LQG quantisation procedure above. However, the 
crucial difference is that the prescription [qHI is tailored to reproduce the standard Fock space 
quantisation in the observable sector, whereas LQG does not. Also, it is not known how to extend 
the prescription of 'o^] to full QED or Yang- Mills theory, or how to implement the procedure in 
a perturbative approach (with unphysical intermediate states in Feynman diagrams). 



El 



^ r 



< CXD 



(4.8) 
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4.2 Comparison with the loop approach in gauge theories 

The (Wilson) loop approach has been advocated in the past in the context of gauge theories, 
as an attempt to non-perturbative quantisation. This approach eventually did not lead to the 
expected success, see jSS, ^] for detailed reviews of the subject. Nevertheless LQG practitioners 
argue that their approach, despite its similarity to the loop approach in gauge theories, should 
work for gravity. Because this relates to one of the frequently asked questions about the LQG 
programme, we here briefly contrast some key features of the Wilson loops in gauge theory with 
the behaviour of spin networks in LQG. Not unexpectedly, the essential difference turns out to 
lie in the special scalar product (|4.7j) used in LQG. 

The first difference between Wilson loops and spin networks is that, due to the product ()4.7() . 
any spin network is orthogonal to any other spin network except itself, including the trivial empty 
one {i.e. the vacuum). Hence, 

{^r,{c}[A]) = {l\^r,{cy[A]) = 0. (4.10) 

By contrast, in gauge theory, the expectation value of a single Wilson loop encodes key physical 
information: it determines the behaviour and properties of quark - anti-quark potentials, and can 
signal confinement of quarks. More specifically, in a perturbative approach we have 



(Wc) = { exp 



ie j) Afj^dx'^ ^ = exp —^j)j)dx^dx''An,y{x — y) 



(4.11) 



where Af^i, is the gluon propagator. In conventional quantum field theory (QCD), this expression 
is divergent, in contrast to ()4.1U() ) because of the self-interaction of the quarks and gluons prop- 
agating in the loop. The way in which it diverges generically depends on the shape of the loop. 
Nevertheless, the UV divergences in ()4.11() can be consistently removed, for example by point 
splitting regularisation, and a well defined answer obtained. All loops can be regularised in this 
way (while it is only known how to renormalise, smooth and non- self- intersecting loops joil). A 
better, and entirely non-perturbative scheme to calculate the W ilson loops and other quantities 
of interest is, of course, provided by lattice gauge theory [ii, il, H • 

In the case of correlation functions of several Wilson loops, no new types of divergences 
appear. For two non-intersecting loops C and C 'self energy divergences' can be removed by 
defining the 'connected correlator' 

{WcWc')conn := {WcWc) - {Wc){Wc') = finite , (4.12) 

which can be thought of as the analog of the scalar product of two single loops in LQG. In first 
approximation we then obtain 



le 
' 2" 



I dx'^ <p dy" Af,^{x - y) 

JCi JC2 



(4.13) 



The main point about this result is that it depends continuously on the shape of the loops Ci and 
C2 - unlike the scalar product or correlator of two spin networks ()4.7() . which is non-zero if and 
only if the two networks completely overlap. Indeed, different shapes of Wilson loops in gauge 
theory are inequivalent, such that the relative change can be expressed via the Yang-Mills fields 
strength and the 'non-abelian Stokes' theorem'. The value of the Wilson loop is invariant under 
continuous deformations only for vanishing field strength (as in a topological theory, such as 
Chern-Simons gauge theory). In LQG on the other hand, the physical states (for which we have 
to enlarge the space S, as we said) are supposed to be diffeomorphism invariant. Hence there is 
no physical information whatsoever in the shape of the loop, since two networks of different shape 
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Figure 4: A flux operator intersecting a three-valent node. The action is given by equa- 
tion 

but identical topology can always be related by a suitable diffeomorphism. This equivalence holds 
independently of the 'value' of the Ashtekar field strength in this state. 

When matter is included in the LQG formalism, consistency would seem to require that 
standard gauge fields must be treated analogously to the Ashtekar connection, i.e. with holo- 
nomies associated to edges, and a scalar product similar to (|4.7|) . leading to a non-separable 
kinematical Hilbert space for the matter sector as well. It is not known how (and whether) this 
formalism can recover known results such as the ones sketched above, and in particular the shape 
dependence of the Wilson loops, and we do not know of any attempts in this direction. At the 
very least, this will require a Sophisticated analog of the 'shadow states' introduced in M- For 
kinematical states (which do not satisfy all the c onstraints) progress along these lines has been 
made recently in pj, |2fl|, 121,, 22,, E3l • 

4.3 The action of elementary operators on a spin network 

We would now like to implement the canonical variables as operators on spin network wave 
functions in such a way that the canonical commutation relations (|4.1j) and (|4.2j) are satisfied. 
This will provide us with the 'building blocks' that will allow us to define the action of the basic 
kinematical operators (area and volume operators) and the Hamiltonian constraint operator on 
any spin network wave function. The elementary holonomy and flux operators defined in (|3.35|1 
and 1)3. 41() act in a rather simple way on the spin network states. Namely, the holonomy {he[A])af3 
simply acts as a multiplication operator, and thus adds the edge e to an existing network F. 
This edge may be disconnected from or 'overlaid' on F; usually it will appear in some gauge 
invariant combination, whose action on the spin network wave function produces another such 
wave function. The action of the flux Fs[E,f] on a network state is either zero (if the surface 
element S does not intersect the graph F anywhere), or otherwise cuts in two any edge which 
pierces S, and inserts a new intertwiner vertex Tq into the spin network at the point of intersection. 
Consequently, on such a single edge (which is part of some spin network wave function \I'r,{c}) 
one has 

F§[E] ( . . . {hM])^^ ■■■)= 8vri/|>% i{e, 5) ( . . . {h,, [A] [A])^^ . . . ) . (4.14) 

where e = ei U 62 and we have now omitted the symbol pj^. Dots indicate the remaining part of 
the spin network wave function. 

If -^^[-E'] cuts the network at a vertex, the result depends on the choice for the intertwiner as 
well as the position and orientation of the surface S with respect to the edges. Let us illustrate 
this by considering the 3-valent vertex depicted in figure |3 The flux acts at the end of every 
edge that emanates from the node. Edges which are located above the surface S contribute with 
opposite sign from the ones that are located below, and for the example of the 3-valent vertex 
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the net result is 



X ((^e,[A])^^^^(/.,,[A])^^^^(/.,3[^])^^^J (4.15) 

(the uncontracted indices (3i, connect to the remainder of the spin network wave function) . 

In the last step, we have used the invariance property of the Clebsch-Gordan coefficients, 

Likewise, when the node has higher valence, the action of the flux operator Fg[E\ inserts 
factors ±pjg(r") in the appropriate place in the spin network wave function, with a '+' for edges 
e 'above' 5, and a '— ' for the edges 'below' S. This action is sometimes symbolically represented in 
terms of angular momentum operators J^" . Consequently, the flux operator can be re-expressed 
as a sum of angular momentum operators, viz. 

e 

where the sum runs over all edges entering the vertex v, and the signs cr(e, 5) = ±1 are determined 
by the position of e in relation to 5. When the signs are all '+', the sum vanishes for gauge 
invariant wave functionals by angular momentum conservation, i.e. Gauss' law. 



4.4 The kinematical operators: area and volume operators 

Two important kinematical operators in quantum gravity are the volume and the area operators. 
The area operator is relevant in the application of LQG to black hole physics, while the volume 
operator is essential for the very definition and the implementation of the scalar (Hamiltonian) 
constraint. It is important to stress that neither of these operators is an observable a la Dirac, 
in the sense that neither commutes with all the canonical constraints (classically or quantum 
mechanically) ^. Their supposed physical relevance is based on the general belief that proper Dirac 
observables corresponding to area and volume should exist, provided suitable matter degrees of 
freedom are included ^103, 104]. The reason is that, for an operational definition of area and 
volume, one needs 'measuring rod fields' in the same way that the operational definition of 
'time' in quantum gravity requires a 'clock field', and that the true observables are appropriate 
combinations of gravitational and matter fields (which then would commute with the constraints). 
Althou gh t here has been much discussion in the literature about these 'relational' aspects (see 

I30I. l3ll. li^ l and references therein, and [lOSl . Il06| | for a more recent discussion in the context 
of LQG), these issues are by no means conclusive lv se ttled, and we are not aware of any fully 
concrete realisation of this intuitive idea (see also |in7i |). 

Setting aside such conceptual worries, the main strategy in constructing the area and volume 
operators is to first write the classical expressions for these quantities using the basic variables 

®As pointed out in [T^. such observables might perhaps better be called 'perennials', to distinguish them from 
the kinematical 'observables' which commute only with the spatial diffeomorphism generators. 
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(holonomies and electric field), and then promote those into operators. To do this one first has 
to regularise the classical expressions by a discretisation of the integrals |l08l . Il09| | . It will turn 
out that the ambiguity in the way one does this is irrelevant for the area operator, while it does 
affect the result for the volume operator. 

We should point out that the area and volume operators appear to be of no particular relevance 
in conventional geometrodynamics. There, they would simply be represented by multipli catio n 
operators (defining a 'length operator' would also be straightforward, as it is in LQG jlinl |). 
As long as these act on smooth wave functionals, there would be no problem, and issues of 
regularisation and renormalisation, like the ones we are about to discuss, would not even arise. 
On the other hand, on more singular 'distributional' functionals, these operators would not be 
well defined, as they depend non-linearly on the basic metric variables. Whether it is possible or 
not to properly define these operators in geometrodynamics therefore depends on what the 'good' 
and the 'bad' wave functionals are. But this question is impossible to answer in the absence of a 
suitable scalar product and a Hilbert space of wave functionals (see e.g. ji^] for a discussion of 
these and other difficulties in geometrodynamics). 



4.4.1 The area operator 

A two-dimensional surface 5 C S possesses a unit normal vector field N and its classical area is 
given by the determinant of the induced metric hrnn ~ Qmn — NmNn- To compute the area, it is 
however simpler to work directly with the area element (|H.H8j) (cf. subsection \',^A^ . so that 

As[g] = / VdF" • dF<^ . (4.18) 



To convert this into an expression that can act on the spin network wave function, we subdivide 
the surface into infinitesimally small surfaces Sj (with / = 1,...,A^), such that S = U/S*/. 
Next, we approximate the area by a Riemann sum (which, of course, converges for well-behaved 
surfaces S), using 

I VdF» • dF'' w ^JF§^ [E] F^^ [E] . (4.19) 

Therefore 

TV 



As [E'^] = lim J F-^ [E] F-^ [E] . (4.20) 

1=1 

The Riemann sum on the right-hand side of 1)4. 2U() can now, for fixed A^, be promoted un- 
ambiguously to an operator. Namely, when acting on an arbitrary spin network wave function, 
the area operator will receive a contribution each time an edge of the spin network pierces an 
elementary surface Si in the sum (|4.20j) . More precisely, one can use equation (|4.14|) to derive 
(always in the appropriate representations pj^, and assuming L{e,S) = ±1) 

F^^[E]F^^[E] = i\8nllhjf {h,,[A] Ta r'^ K,[A]) 

(4.21) 

= i8nllh-ff {jeUe + l)) {he[A])^^. 

Note also that if the the surface Sj had been pierced by two edges (carrying representations ji 
and j2) instead of just one, the previous expression would generalise to 

F^E]Fl[E] {h,M])^^{heAA])^s 

= {8nllh^f (ie,(ie, + 1) + Je.(je. + 1)) (/^e, [A]) (/j^, [A]) . (4.22) 
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Hence if one applies the operator (|4.2U|) to wave function associated with a fixed graph T with 
L(r) edges, and refines it in such a way that each elementary surface Sj is pierced by only one 
edge of the network, one obtains 




(4.23) 



p=i 



It is obvious that further refinement of the area operator then does not change the result. It is 
also clear that the final result, therefore, does not depend on the shape of the elementary area 
cells. From ()4.23|) . one sees that all spin networks are eigenstates of the area operator (when 
the Si intersect only edges, see below). Due to the 'discreteness' of the spin network states, 
the spectrum of the area operator is also discrete. The quantisation scale is set not only by the 
Planck length Ip, but also depends on the parameter 7.^'^ 

We have so far assumed that the surface elements Sj in the above computation 'meet' only 
the edges (hence insert a bivalent vertex into the spin network), but not the vertices of the spin 
network. Making use of the general result (|4.17)) we can also evaluate the area operator when 
the surface elements Si intersect a vertex. In this case, the result depends on the particular 
intertwiner at that node and on the orientation of the surface element. Let us illustrate this on 

a four-valent vertex, with intertwiner C''^'^^^r,C^^''^''n (see figure 0. Consider first the situation in 

■ 0:10:2/5 0304P ^ 

which the edges 1 and 2 are both located on the same side of the surface element and the same is 
true for edges 3 and 4 (figure EK)- In this case, one can use the invariance of the Clebsch-Gordan 
coefficients (|4.16|) to move around the inserted r-matrices, and re- write the action of the area 
operator in terms of a 'virtual' edge. For our example, this yields the result 



Observe the relative factor of \/2 with respect to (|4.23jl : the result for a virtual edge is not the 
same as the one for a real edge. The four-valent node is thus an eigenstate of this area operator. 
Next, consider the situation in which the surface element is oriented as in figure EJa. In this case, 
one first has to rewrite the product of Clebsch-Gordan coefficients by means of the recoupling 
relation (|4.6|) . Independently of the normalisation of the 6j-symbol, the important thing to note is 
that the area operator will change the relative coefficients of the various terms in the sum ()4.6|) . 
picking up a contribution oc m{m + 1) for each admissible value of the angular momentum 
Iji ~ isl ^ < ji + j3 in the new channel 1 + 3 — > 2 + 4. Hence, the four-valent vertex is 
generally not an eigenstate of the area operator associated to the surface of figure Eb. For higher 
valences, the number of different contributions increases rapidly with the number of accessible 
channels. Therefore, the area operator in general does not act diagonally on an arbitrary spin 
network state. 

We conclude this subsection with three remarks. Although infinite spin networks have been 
excluded hy fiat, it is nonetheless instructive to see what would happen if we did include them. 
In this case it may not be possible to achieve the refinement of one edge per one cell in a finite 
number of steps. To see this, consider the example of a strongly fluctuating network, of 'shape' 
z = y = sin(l/x) in the fiducial background coordinates. It is clear that if one considers the 
area of the y = surface, then no matter how small the size of cells, the cell around zero will 
always be pierced by an infinite number of edges, and hence the area operator would not be 
well defined. Because classical space(-time?) presumably emerges in the limit of infinite spin 

^"We emphasise that this discreteness, which is often invoked as the underlying reason for the absence of 
divergences in LQG, hinges on the compactness of the group SU(2), and would not hold if this group were replaced 
by a non-compact or complex group (as would be the case for the original Ashtekar variables with imaginary 7). 




(4.24) 
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Figure 5: The action of the area operator on a node with intertwiner C^ai ^2/3^ 0304 fs- 
figure on the left, the location of the edges with respect to the surface is such that the invariance 
of the Clebsch-Gordan coefhcients (|4.16|l can be used to evaluate the action of the area operator. 
The result can be written in terms of a "virtual" edge. In the figure on the right, however, the 
edges are located in a different way with respect to the surface. The invariance property (|4.16|) 
does not apply, one has to use the recoupling relation H4.6|l . and the spin network state is 
therefore not an eigenstate of the corresponding area operator. 



networks, the exclusion of such 'ih-behaved' networks reaUy amounts to a prescription for the 
order in which limits are to be taken, when they would otherwise produce an ambiguous answer. 

Secondly, the area operator is unbounded, and not every state in "Hkin possesses finite area 
(whereas every element of S has finite area) . This can happen even for compact spatial manifold 
E, and the wave function can then be thought of as associated to an 'infinitely crumpled surface' 
in E. To see this more explicitly, consider a sequence of spin network states with edges 
labelled by p = 1, . . . , L(r„), and spins jp"^; then we have 




K^^E I E + l)a„^„ I . (4.25) 



=1 \ p=i 



The expectation value of As can be made to diverge by letting the spins {jp^^} grow sufficiently 
fast with n, even if |a„p||^'n|P < 00. 

Thirdly, the formula 1)4. 23() plays an irnportant role in the microscopic explanation of the 
Bekenstein-Hawking entropy of black holes 0, |^ , which is now reduced to a counting problem. 
As an important prerequisite for this calculation, one must first define the notion of (isolated) 
horizon in the quantised theory, i.e. introduce a condition on the quantum states that singles out 



the 'horizon states'. This has been accomplished in |43]. Fitting the solution of this problem to 



the known factor relating entropy and area of the horizon then yields a prediction for the value 
of the parameter 7. While the original value given was 7 = 0.12738402 . . ., a more recent 

calculation gives 7 = 0.23753295 . . . [ll3,[ll2 (see also Elffll). 



4.4.2 The volume operator 

The construction of the volume operator is less 'clean cut' than that of the area operator, be- 
cause it is fraught wit h am biguities, which can only be eliminated by making certain choices 
and by averaging, see for a comprehensive discussion. For a different der ivati on using a 

more conventional point-splitting method, which yields the same final result, see [ll5l|. In a first 
approximation one might try to use the functional differential operator (|3.3H) . but its direct appli- 
cation would lead to (square roots of) singular factors (J*-^-* (x, x). The 'detour' in the construction 
to be described, in particular the smearing over the surfaces Sf in (|4.30j) below, is necessary 
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Figure 6: An elementary cell in the Ricmann sum (|4.29|l . together with the three surfaces on 
which the fluxes are evaluated. 



mainly in order to eliminate these singular factors. One starts with the classical expression for 
the volume of a three-dimensional region C S, 



v{n) 



(4.26) 



Just as with the area operator, one partitions O into small cells = Uji^i, so that the integral 
can be replaced with a Riemann sum. To express the volume operator in terms of canonical 
quantities, we rewrite the volume element dV — 
dF^ = 6*2 A d^,. . . (see (lOHll l as 



6 A 9 A 6 in terms of the area elements 



dV = ^/\dF^^dF^^dF^ 



(4.27) 



In order to express this volume element in terms of the canonical quantities introduced before, 
we next approximate the area elements dF°' by the small but finite area operators such 
that the volume is obtained as the limit of a Riemann sum 



lim 




1=1 



(4.28) 



Now, in order to make sense of this expression, we must for each cell $7/ choose three small 
non-coincident surfaces Sf, which subdivide O/, as shown in figure IHl This should be done in 
such a way that the r.h.s. of (|4.28j) reproduces the correct classical value. For instance, one can 
choose a point inside each cube fi/, then connect these points by lines and 'fill in' the faces. In 
each cell 0/ we then have three lines labelled by a = 1, 2, 3; the surface Sf is then the one that 
is traversed by the a-th line. With this choice it can be shown that the result is insensitive to 
small 'wigglings' of the surfaces, hence independent of the shape of Sf, and the above expression 
converges to the desired result 

Having constructed a sequence of refinements with the correct classical limit, we can now 
attempt to promote (|4.28|) to a quantum operator by invoking the known action of each factor Fga 
on a spin network wave function, viz. 



N 



v{n) 



lim 



N- 



lim 



7=1 



N- 



N 

,E 

7=1 



qf\S},SlSf)\ 



(4.29) 



Alternatively, one can partition E with other types of cells, with analogous results. 
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Figure 7: A series of refining steps for one cell in a Ricmann sum, such that the indicated node 
is never at the centre of any of the boxes except in the strict N ~^ oo limit. A volume operator 
based on cells of this kind would always yield a vanishing result for any finite cell size. 

with the operator (= square of the local volume) 

qf\s},SlS]) := ^ea,c^[E]F's][^]%[E]. (4.30) 

Ignoring the square root for the moment, (|4.3n|) is indeed well defined as an operator at this 
point. However, unlike the classical volume, the action of q^j^^ on the spin network wave function 
depends on the way we partition Q, and more specifically on the position of the faces Sf in relation 
to the edges and vertices. As a simple example of the inherent ambiguities consider a cube which 
is traversed by one edge e, but not containing any vertex of the network. If the surfaces Sf are 
positioned in such a way that the edge intersects each one of them with t(e, Sf) 7^ at three 
different points in Qj, the operator simply inserts the matrices T^,r^ and at the appropriate 
places in /ie[^]- If the cell is sufficiently small, we can bring these matrices together, using 
and thus pick up a contribution of order one. It is easy to construct sequences of 
refinements such that an arbitrarily large number of cells fi/ gives such a contribution, no matter 
how small the cells are. In order to prevent catastrophic divergent sums from appearing on the 
r.h.s. of H4.29() one must therefore rule out such configurations 'by hand' and stipulate that each 
cell 0/ may contain at most two such intersections with any edge, in which case one of the three 
factors in (|4.3U|) will give zero. 

With this proviso, the sum receives contributions only from those cells which contain 
a vertex v. If the vertex v does not coincide with the intersection point of the three faces Sf, 
further refinement will yield more cells, which do not contribute by the above requirement. Hence, 
with sufficient refinement only those cells will remain and contribute for which the vertex v sits 
precisely at the intersection: v = SjCiS^nSj. This means that we must also exclude refinements 
such as in figure [7| for which SjnSjn Sj 'misses' v by an amount that gets smaller with N, but 
never vanishes for N < 00. With these assumptions, the resulting action of the operator qj on 
the vertex v reduces to 

qj = {87rlhfj^ e{eue2,es)eatJ!,'^Ut'Ut\ (4.31) 

61,62,63 

where the sum runs over all possible triplets of edges in the vertex v, and the angular momentum 
operator (defined in H4.17() l by definition inserts a matrix Pj^iT"") at the vertex v into the 

spin network wave function. The coefficients 6(61,62,63) are orientation factors which are equal 
to ±1. They depend on how the faces Sf cut across the vertex, and on which edges are 'above' 
and 'below' Sf. More precisely, in terms of the factors a{e, S) introduced in (|4.17j) they are given 
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Figure 8: Two different choices for the relative orientation of the edges emanating from a 
four-valent node and the cell used in the definition of the volume operator (prior to averaging 
over the positions of the surfaces Si). It is always possible to choose the division of the cell such 
that all four edges emanate from the same octant. We have here drawn the 'active' version in 
which the edges rather than the surfaces are deformed. Note that these two configurations are 
not homotopically related. 

by 

e(ei, 62, es) = a{ei, Sj) a{e2, S|) a(e3, Sj) (4.32) 

and thus completely fixed by the relative positions of the edges and the surfaces. Let us also 
stress they are independent of any fiducial metric background structure, and purely combinatorial 
- hence invariant under diffeomorphisms. ()4.32|) is not yet the final answer, but let us first make 
the above formulas a little more concrete by working out the action of the above operator on a 
4-valent vertex (|4.5|) . Its net effect is the replacement of the intertwiner according to 

CQV.V.a4 ^ const. ^ 6(61,62,63) eabc 

perm.{l,2,3,4} 

" («-.<^°>)..,. («-<^'>)„,, («-<^'>)„* ^"•'"^S':^ . 

where the sum runs over all choices of three edges out of four. From this formula it is obvious 
that in general the operator qj acts non- diagonally on a spin network wave function |r, {j}, {C}), 
because it changes the intertwiners {C}, although it does not affect the spin network itself {i.e. the 
graph r). The operator vanishes when all the edges enter v through the same octant, in which 
case e = +1 for all choices, and 

eabc{JtJlH + JUlJl + JUPl + J2JIJD = (4.34) 

by angular momentum conservation, Jf + J2 + Jf + = 0. Note that we can in principle arrange 
any 4-valent vertex to conform with this choice (see figure IB} • Similar results hold true when yet 
more edges enter the vertex from one octant. 

Let us now return to the coefficients £(61,62,63) in (|4.32j) . Though manifestly background 
independent they do depend on how we choose to position the faces and edges: for instance, 
there is no a priori argument that would force a given edge linking two vertices in adjacent cells 
to appear 'to the left' or 'to the right' of the face connecting the two vertices, as illustrated 
in figure IHl and in fact almost any combination of factors ±1 is possible. To eliminate this 
remaining ambiguity in the definition of the volume operator, one now averages over the various 
possibilities More specifically, this is done by integrating over the positions of the surfaces 

^^We might note that this averaging is 'undemocratic' insofar as it is applied only to cells containing a vertex. 
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Figure 9: Sign ambiguities in the volume operator, present because the relative location of the 
three surfaces Sj and the edges of the spin network can be chosen in different ways (again prior 
to averaging, as in figure |HJ|. 

with some measure d/i(0i, 02, ^3), where the 9i are suitable angular coordinates, viz. 

e(ei,e2,e3) := J d^(0i,02,03) <^{ei, S}{9i)) a{e2, S]{e2)) <T{e3, SjiOs)) . (4.35) 

This leads to a well defined combinatorial answer depending on the relative orientation of the 
three edges (ei, 62, 63) (= the sign of the determinant of the dreibein spanned by the tangents to 
the three edges at the vertex) , but only modulo an overall factor which encapsulates the freedom 
in the choice of the measure d/i. In principle such a factor can arise for each vertex, but in order 
to preserve diffeomorphism invariance of the volume these factors should be chosen to be the 
same for all vertices. This factor is the one ambiguity w hich remains in the definition of the 
volume operator, and this ambiguity cannot be eliminated |l09l |. 

Finally, given an assignment of e{ei,e2,es) for each vertex, how do we actually compute the 
spectrum of the volume operator in (|4.29l) ? First of all, it is easy to see that this operator always 
vanishes on three- valent vertices which satisfy the Gauss constraint. This is a direct consequence 
of angular momentum conservation, or equivalently 

(^°')„.. (-'')„., (ft... (-'')„.. ^^^t' = » ■ ("6) 

In other words, the volume operator vanishes on gauge invariant two- and three- valent vertices 
(but not on gauge variant ones), and picks up contributions only from spin networks with vertices 
of valence > 4. Given any spin network wave function, the first step is to work out the action 
of qi according to the above rules. In a second step this action must be diagonalised for each I on 
the spin network wave functions. This can be done because the operators qj commute amongst 
each other, i.e. [qi,qi'] = for / / I'. Only now can we define the absolute value \qi\ and 
take the square root -^/[^/T using the spectral decomposition of this operator (note that qj is not 
necessarily a p ositiv e operator). Unfortunately, these steps are all but straightforwar d to carry 
out in practice |lld |. The analysis can be simplified by applying the operator identity lisl . 117 1 



eatcJlJir^ = \ [(J^ + J^)^ (J^ + J^')^] (4.37) 

to any three edges entering the vertex, enabling a numerical analysis of simple configurations with 
low- valence nodes. This procedure has recently been used to show that the operator corresponding 



If it were extended to other cells, some of the unwanted features which we just disposed of might re-appear. 
Vertex-less cells which satisfy the condition that edges intersect at most two surfaces do not necessarily satisfy this 
condition for all surfaces orientations which appear in the averaging procedure. 
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to the inverse volume is not bounded from above when acting on zero-volume spin network 
states ji^l (in marked contrast to the symmetry reduced models of [s^, where the avoidance of 
an initial singularity hinges on the fact that the operator corresponding to the 'inverse volume' 
is bounded from above^^). Nevertheless, the resulting formulas remain very complicated. The 
lack of a more explicit representation for the operators \fqi is also a serious impediment towards 
a better understanding of the Hamiltonian constraint. 

4.5 Some comments on regularisation dependence 

At this point some comments regarding the regularisation procedures for kinematical operators 
which we have outlined above, are in order. Similar comments will apply to the regularisation 
of the Hamiltonian constraint operator. First of all it should be evident from the foregoing 
discussion that the notions of 'finiteness' and 'regulator independence' as currently used in LQG 
are not the same as in conventional quantum field theory. There one insists that renormalised 
quantities should come out the same (modulo a finite number of normalisations) after removal 
of the regulator, no matter how the latter is chosen. By contrast, LQG does not claim that 
different prescriptions for constructing an operator will always lead to the same physical result. 
This is well illustrated by the volume operator, whose definition and whose very existence hinge 
on various choices that have been made along the way. 

Let us nevertheless briefly summarise what is meant here by 'regularisation'. In order to 
regulate a given classical phase space quantity defined by an integral (such as the volume), one 
first approximates that integral by a Riemann sum (which is of course assumed in the classical 
theory to converge to the desired quantity in the limit of infinite refinement). One then replaces 
in this finite Riemann sum all terms by operators, usually in such a way that the various area 
or volume elements are absorbed into the definition of these operators (as exemplified by the 
'small' operators Fga[E] entering the definition of the area and volume operators). Acting with a 
sequence of these discretised operators on an arbitrary, but fixed spin network state, the resulting 
sequence of refinements (labelled by N, say) eventually becomes stationary (i.e. constant for N 
sufficiently large), because the state is kept fixed as the limit — > oo is taken. This follows 
straightforwardly for the area operator, but requires a choice of 'preferred' regularisations for the 
volume operator. With these choices, the 'correct' result is always obtained already after a finite 
number of refinements. A good analogy would be the Riemann integral of a step function, for 
which a suitably refined Riemann sum attains the exact value already after finitely many steps, 
after which the result no longer depends on the 'regulator' (i.e. the size of the intervals defining 
the Riemann sum). In this sense, the results of LQG are indeed 'regulator independent'. 

Prom what has just been said, it is evident that infinities can never appear in the LQG 
regularisation procedure, and in this sense the resulting theory is 'finite', at least as far as the 
kinematical operators are concerned. LQG nevertheless requires the regularisation of the area 
and volume operators in order even to be able to define the quantum counterparts of the classical 
constraints via Thiemann's trick. The outlined regularisation is, therefore, not introduced to 
remove divergences. Standard short distance, QFT divergences 'disappear' in the LQG approach 
by the very construction of the theory: all states are discrete, and at any step of the calculation 
one deals only with a finite number of objects. The price one pays are ambiguities of the type 
encountered above, some of which can only be eliminated by making ad hoc choices, such as the 
rule prescribing the way in which an edge has to traverse a cell. We shall encounter more such 
ambiguities when we attempt to define the Hamiltonian constraint operator. 



^■^The operator X corresponding to the inverse volume is not the same as the inverse of the volume operator. 
The latter does not exist on all of Tikin because of the occurrence of zero eigenvalues in its (discrete) spectrum. 
Instead, the operator X obeys XoV = l + 0{h), and can therefore in principle stay bounded 'near the singularity'. 
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Of course, these considerations do not necessarily imply any inconsistency. Namely, one can 
always adopt the above prescriptions for the kinematic operators as definitions, which can be 
plausibly related to the corresponding classical objects, and in this sense the definitions of the 
above operators in LQG are perfectly acceptable. What we would like to emphasise, however, 
is that with this construction the work is hardly begun. The far more difficult task is now to 
show that these definitions, when put together and inserted in the constraints, yield a sensible 
quantum theory, and a consistent quantum constraint algebra, in particular. The situation here 
is reminiscent of the one in the early days of constructive quantum field theory (see e.g. |lis| | 
and references therein). The problem then was not so much finding non-perturbative regulators 
and establishing the existence of selected limits when the regulator is removed. Rather, the main 
problem was to establish that the resulting theory has all the desired physical properties, to wit, 
full Poincare invariance, locality and causality. As is well known, no interacting quantum field 
theory in four dimensions satisfying all these requirements is known to this day And for 
those low dimensional models which could actually be shown to exist and possess all the requisite 
properties, the tightly knit construction left no room for ambiguities other than the normalisation 
of a finite number of physical parameters (masses and coupling constants) . 



4.6 Coupling to matter fields 

Any realistic theory of quantum gravity should at least allow for a coupling to the matter fields 
of the Standard Model of elementary particles. In the framework of geometrodynamics it is fairly 
easy to extend the relevant formulas to incorporate matter fields (but, of course, the difficulties 
with quantisation remain the same as before). Likewise, it is possible to include gauge fields and 
fermions in the LQG formalism, essentially without any restriction on the gauge group or the 
number and structure of fermion families — as we already pointed out LQG does not appear 
to put any restrictions on the matter couplings (such as renormalisability) . Armed with the 
translation rules of section [3.21 the transition from the 'old' metric or dreibein variables to the 
loop variables is in principle straightforward, though usually cumbersome (and, we may add, of 
considerably diminished esthetical appeal in comparison with Ashtekar's ori ginal reformulation 
of pure g ravitv). Because there is an ample literature on this subject (see e.g. 



[l2J]), and because the introduction of matter couplings does not alleviate the problems with 
the Hamiltonian constraint in LQG, we will be very brief and only mention two selected aspects 
here. 

Consider, for example, the coupling of fermion fields in the LQG setting [llfll . Il25l |. The 
starting point is the classical continuum action of a Dirac fermion coupled to gravity 

S = jd^Exi^EA^V^x^ jd^xExi^EA^{d^ + u:^BCl'''')x. (4.38) 

where we have omitted a possible mass term for simplicity, as its inclusion would pose no new 
conceptual problems. There are two immediate questions one faces when trying to construct 
quantum Hamiltonian for this system. One is the choice of the conjugate pair of variables for 
fermions; finding the proper fermionic variables may also require redefinitions of the other fields 
so as to 'diag onal ise' the canonical (Poisson or Dirac) brackets. For instance, for the above model 
it turns out |l2(]| | that a good choice for the conjugate pair of fermionic variables is the pair of 
half-densitised spinors 

X = V det E il) , 7r = ixh°. (4.39) 



^ Even the best candidate model for realising the quantum field theorist's dream of a UV finite theory, A'^ = 4 
supersymmetric Yang-Mills theory, has not been rigorously shown to exist beyond perturbation theory. 
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If instead of densitised spinors we used bare spinors (without any facto rs of the determinant of 
the dreibein), one would have to modify the gravitational connection jl2Clj | in order to obtain 
diagonal canonical brackets, but at the expense of making the connection complex again 

A'in = - he + iK^^ + ^e^xx ■ (4.40) 

^ /I A / J? 



This leads to a new reality problem, similar to the one already encountered in 76], and not 
different from the one which eventually led LQG to abandon the 'old' Ashtekar variables with 
7 = ibz. The second question is how to express all variables in terms of the LQG variables 
he [A\ and [E] . The answer to this question can be obtained by employing the identities (|3.26|1 
and H3.27() . A similar procedure will be spelled out in the next section in order to construct the 
Hamiltonian constraint, so we will refrain here from giving more details. 

When implementing the matter Hamiltonians in terms of spin networks, the description of 
(abelian or non-abelian) gauge fields and matter fields is somewhat similar to the description of 
standard lattice gauge theories: matter fields (scalars or fermions) are attached to the vertices of 
the spin network, while gauge fields ar e ass ociated to the edges (links), see e.g. (121. . .122 . .123] for 
a treatment of electromagnetism and |l24l | for scalar field theories in the LQG formalism. There 
are, however, two main differences between lattice gauge theory and LQG. In the former, the 
space or spacetime lattice is given, and hence there are only three (or four) possible directions for 
the edges. The coupling between two fermions xi and X2 on adjacent vertices 1 and 2 connected 
by an edge e = (12), is uniquely given by Xi7i2^i2[^]X2! where 712 is the 7-matrix 'along' the 
edge (12). In LQG (and in general relativity), this simple formula fails because the direction 
of the link is arbitrary, and a vielbein factor must be inserted to 'align' the 7-matrix with the 
holonomy. When rewriting the dreibein by means of the formulas of section 13.21 this usually 
leads to rather awkward expressions in terms of the spin network variables he\A\ and Fg[E\. 
This problem is further compounded by a second difficulty which goes back to the issue of what 
the correct value is for the Barbero-Immirzi parameter. As we mentioned already, for 7 7^ iti, 
the connection A is not the pullback of the (chiral) spin connection. Therefore, the holonomy 
/ii2[^] does not effect the correct parallel transport between xi and X2, and must accordingly be 
amended by further (and again rather awkward) modifications. 

We should stress that, even when one has obtained an expression for the matter action in 
terms of spin network operators, one does not know whether the resulting quantum theory has any 
relation to the usual Fock-space quantised theory. Showing that such a relation exist requires the 
field-theory analogues of the 'shadow states' for point particles These states are complicated 
linear combinations of elementary spin network states, approximating Fock s pace coherent states. 
In flat space some results have been obtained in the case of Maxwell Il2ll. Il22l. Il2,'l and scalar 



fiel ds 11241 . Rece nt at t emp ts to construct coherent states in the presence of gravity include [1 
20,|2l|,y,|2i,|2i,[lo],|lo3. 



A completely different question concerns the role of the structure group SU(2), and the 
possibility of generalising Ashtekar's variables to higher dimensions, to achieve a Kaluza-Klein 
type unification of matter and gravity. The closure of the classical constraint algebra in the 
Ashtekar reformulation of gravity requires use of the identity f-ahe^^^'^'^ = 25^^, which is only 
valid for the structure constants of the group SU(2) and its non-compact form S0(l,2). These 
properties no longer hold for SU(A^) when N > 2. For this reason, the LQG reformulation of 
gravity works only in thr ee an d four spacetime dimensions. Nevertheless, one can try to generalise 
the formalism to SU(iy) 112611 and thereby arrive at a novel type of coupling between gravity and 
Yang- Mills fields 1271 . Il28l |. A very different attempt to generalise the formalism to include 
matter couplings is based on the observation that an inve rse d ensitised vielbein also appears 
naturally in certain reformulations of I? = 11 supergravity [l29l |. It is not known whether any 



34 



of these ansdtze can be implemented in a spin network formulation (although replacing SU(2) 
holonomies by SU(A^) holonomies seems straightforward enough). 

5 Quantum constraints 

As we have seen in section [3.11 the Hamiltonian framework leads to three classical constraints: 
the Gauss, diffeomorphism and Hamiltonian constraints. To implement these constraints at the 
quantum level, one must first properly define them, i.e. express them in terms of the elementary 
variables, the holonomies and fluxes, and then investigate their properties. We will first discuss 
the kinematical constraints, which are 'easy' (relatively speaking), whence we can be brief. We 
will then turn to the truly difficult part of the story - the Hamiltonian constraint. 

5.1 Kinematical constraints 

The Gauss constraint is already expressed in canonical variables he[A] and Fs[E, /], and one can 
solve it exactly. This is simply because fulfilment of the Gauss constraint is equivalent to the gauge 
invariance of the spin network wave function: invariance is simply achieved by contracting the 
SU(2) representation indices entering a given vertex in an SU(2) invariant manner. Equivalently, 
one says that angular momentum is conserved at each vertex, viz. 

e 

for each f G F. Recall (see the discussion around ()4.17|) ) that the operator jjf'^^ is defined to 
act by inserting a matrix Pj^ir"") at the vertex v in the representation belonging to the (oriented) 
edge e, with appropriately contracted SU(2) representation indices, and appropriate orientation. 
For the following discussion, we will tacitly assume that the wave function satisfies the Gauss 
constraint, whenever appropriate. 

The diffeomorphism constraint is more difficult to impose — unlike in geometrodynamics, one 
cannot immediately write down formal expressions which are manifestly diffeomorphism invariant, 
because the spin network functions are not supported on all of S, but only on one-dimensional 
links. Nevertheless, it would seem most natural at first sight to try to solve the diffeomorphism 
constraint on the kinematical Hilbert space Wkin- However, this is impossible for two reasons. 
First of all, the diffeomorphism generator does not even exist as an operator: there is no way 
to differentiate finite diffeomorphisms so as to obtain 'infinitesimal' ones due to the lack of weak 
continuity (and a fortiori, differentiability). Secondly, with the exception of the trivial state 
^' = 1 (the empty spin network), there is just no diffeomorphism invariant state in Wkin- This is 
intuitively clear because every spin network state changes under a diffeomorphism which moves 
the associated graph around. For these reasons, LQG invokes a variant of the so-called 'group 
averaging' method, but in a way that is different from the standard group averaging procedure. 
Namely, application of the latter would require knowing an integration measure d/x(0) on the 
infinite dimensional group Diff(S). The formally diffeomorphism invariant state would then be 

^r°;[^]= / d/i((^)^r,^[Ao</>], (5.2) 
JDifr(E) 

and if the measure existed, this state would legitimately belong to Wkin- ^or the infinite dimen- 
sional group Diff(S) no such measure is known (the 'counting measure' (|5.6j) introduced below 

^^This problem disappears altogether in a topological theory, such as (2+l)-dimensional gravity, where the 
connection is fiat and diffeomorphism invariance is manifest from the outset. 
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does not 'live' on Wkin; but leads out of the kinematical Hilbert space), hence the above expression 
is devoid of meaning — confirming the conclusion above that Wkin does not admit non-trivial 
solutions. LQG circumvents this difficulty by implementing 'group averaging' in a way that is 
more adapted to the non-separable Hilbert space and the properties of the scalar product ()4.7() . 
The key observation here is that the constraint can be formally solved on a much larger space, 
namely the algebraic dual D* of some dense subspace T) of Wkirn such that 

V C Wkin C V* . (5.3) 

Typically, the dual space V* is thus a space of distributions. This is not as outlandish as it may 
seem at first, because for c onstrained systems, the solutions of the constraints are generically 
distributions, see e.g. |l3Cll ] Given any (densely defined and closable) operator O and its 
adjoint with domain T) C Wkinj the action of the dual operator O* on the algebraic dual V* 
of T) is defined for any particular element X G D* hy requiring 

{0*X\^) = {XIO^"^) (5.4) 

to hold for all ^ £ D. The space P* is very large, and in practice one restricts attention to a 
physically motivated vector subspace thereof, the so-called 'habitat' V* G T>* , whose specification 
is closely related to the definition of the Hamiltonian constraint. Observe also that the above 
definition does not require the operator O to be hermitean. 
The obvious ('canonical') choice for D is 

V = S (5.5) 

that is, the space of finite linear combinations of spin network states. However, there may exist 
other viable choices for D. The dual space S* is much larger than Wkim and actually too large for 
our purposes. What is important here is that S* contains a certain subspace of diffeomorphism 
invariant states, and its completion w.r.t. to the norm (|5.8|) below. This subspace is contained 
in any of the proposed 'habitats' V*, and is obtained by averaging spin network states. More 
specifically, for any ^'r G S associated with a particular graph T we define the 'group average' 
by the formal sum 

^(^fp) = := (j)* o (5.6) 

<^eDifr(s|r) 

where ((/>* o \E')[A] := '^[A o and where Diff(S|r) consists of the subgroup of Diff(S) obtained 
by dividing out the diffeomorphisms leaving invariant the graph V on which ^'r lives. While 
^'r G 5 by assumption, averaging the state in this way throws it out of S and also "Hkin- At first 
sight, this seems to make matters worse: we have replaced the ill-defined integral over Diff(S) 
by a highly formal continuous sum. However, the day is now saved by two facts: (i) unlike the 
integral, the sum no longer requires unavailable detailed knowledge of the properties of Diff(S), 
and {a) when applying an averaged spin network state ??(4'r') (considered as a bra- vector) to 
any other spin network state ^'r[^], the sum 

{r]{^Y')\^r) = X] (0*o^r'|^r) (5.7) 
</)eDifr(s|r') 



^^A simple example is the solution of the Hamiltonian constraint for a relativistic point particle - the Klein 
Gordon equation - which is — S{p — m ). As is well known, a scalar product and a Hilbert space can be 
defined on this constraint hypersurface, but the scalar product is not positive definite. These features c an also be 
made comple tely explicit in more complicated exactly solvable theories, such as pure gravity |l3ll . Il3^ and pure 
supergravity |l33|| in three dimensions. 
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consists at most of a finite number of terms, and it is this fact which ensures that T/(^'r) is indeed 
weh defined as an element of the dual space S*. Namely, there is no contribution at all if T and 
r' are not diffeomorphic. For T and T' diffeomorphic, only that term contributes for which T and 
0* o r' coincide. In the above equation, the bracket on the r.h.s. is just the scalar product ()4.7() . 
while the bracket on the l.h.s. denotes the dual pairing between S and S*. 

We can now promote the space of averaged diffeomorphism invariant spin network states to 
a pre-Hilbert space by 'dividing out the volume of the gauge group' ( = Diff(S)). More precisely, 
we define the scalar product between two such states as 

iivi^r',^') I V{^r4'))) ■= (^(^r',^) | %,v>> • (5-8) 

Here the bracket on the r.h.s. is the dual pairing between S and S*, and the bracket ((.|.)) on the 
l.h.s. denotes the new scalar product. The new Hilbert space TCdis of diffeomorphism averaged 
states is then obtained by completion w.r.t. this norm, and can be thought of to consist of 
the averaged versions of the states (|4.8|) (but let us emphasise again that it is not a subspace 
of Wkin)- Observe that the Hilbert space Wdiff is still a 'small' subspace of S* because the vast 
majority of elements in S* does not correspond to diffeomorphism averaged spin network states, 
and cannot be obtained as limits of such states w.r.t. to the above norm. Moreover, Tidm is 
'smaller' than the original kinematical Hilbert space TChin, because distinct elements of the latter 
which are related by a diffeomorphism get mapped to the same element of S*; equivalently, the 
averaging map has a very large kernel. There is no way to extend the above norm from Wdiff 
to all of S* , in accordance with the fact that distribution spaces cannot be made into Hilbert 
spaces, but at best into topological vector spaces; however, there appears to be no consensus as 
to which of the many possible top ologies is the right one for S* . While Wkin is non-separable, it 
is expected that Wdiff separable 1341 . 135^ provided one averages over homeomorphisms rather 



than diffeomorphisms (if one divides out only by analytic or C°° diffeomorphisms, TCdm remains 
non-separable, see also for a recent discussion of this point). 

Which, then, is the arena, or 'habitat', in which quantum gravity takes place according to 
LQG, and where one must ultimately consider the action of the Hamiltonian constraint? The 
foregoing discussion shows that there is no easy answer to this question, and apparently also no 
complete consensus in the community as to which space is 'right'. The Hilbert space TCdis of 
diffeomorphism averaged states certainly appears to be a good starting point for implementing 
and solving the Hamiltonian constraint, but it may not be good enough. The reason is that the 
action of the Hamiltonian constraint in general does not preserve 'Hdiffj in other words, it maps 
diffeomorphism invariant states into ones that are not invariant. Consequently, when checking 
the commutator of two Hamiltonian constraints, one must again pass to a larger space, which 
is preserved by the action of the Hamiltonian, but which is no longer a Hilbert space. An 
explicit example are the 'vertex-smooth' states introduced in [2^ . obeying (A'l^'i) = (A'|^2) for 
all ^'i, ^'2 which are related by a diffeomorphism that leaves the vertices of the spin network fixed. 
Alternatively, it has been proposed to implement the Hamiltonian constraint on the kinematical 
Hilbert space Wkin via a weak* limit 0|, see also the following section. At any rate, the challenge 
here — as with any other constrained quantum mechanical system — is to find a subspace that is 
annihilated by all the constraints, and to define a physical inner product which yields the 'final' 
physical Hilbert space. 



5.2 Hamiltonian constraint 

As we explained, the diffeomorphism constraint is imposed not in terms of an operator constraint, 
but through the formal group averaging method. By contrast, the Hamiltonian constraint is 
imposed as an operator constraint. Hence, as a zeroth step, one has to construct this Hamiltonian 
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operator. As we explain in more detail below, in this derivation one is forced, at a certain stage, to 
make choices concerning the action of the Hamiltonian operator on a generic spin network state. 
Thus the resulting definition is indeed motivated by the classical expression for the Hamiltonian 
(expressed in loop variables), but it is plagued by a very large number of ambiguities To test 
the proposed prescription and compare it with other possibilities, one would need to find states 
which are annihilated by the Hamiltonian, and check that in the semi-classical limit these states 
produce sensible results. This is, however, currently out of reach: due to the great complexity 
of the proposed Hamiltonian, not a single physically interpretable eigenstate is known. While 
this should not come as a surprise (after all, we have no reason to expect to be able to solve the 
theory exactly), it is remarkable that even the far more modest goal of working out in complete 
detail the action of the (regulated) Hamiltonian on a general spin network wave function appears 
difficult to attain. 

Even if we can work out the action of the Hamiltonian on any given element of S, we are 
not done. As explained in the previous subsection, in order to ensure satisfaction of the diffeo- 
morphism constraint, we must transfer the action of the Hamiltonian to the dual space S* , or 
at least to some physically motivated subspace V* thereof. Below, we will be mostly dealing 
with a regulated Hamiltonian H[N,€] (where is the lapse function). Here, the regularisation 
parameter e > enters via a plaquette -P(e), which is attached to a vertex and which must be 
shrunk to zero at the end of the calculation. However, for a given ^, the limit e — > of H[N, e]^ 
does not exist on S, because wave functionals supported on the same network, but with an extra 
loop dP{e) attached to one of the vertices, are orthogonal to one another for different values of e 
by 1)4. 7|) . For this reason, one must resort to a weaker notion of limit by transferring the action 
of the Hamiltonian to the dual space. More specifically, we define the limit e by demanding 



to hold for all "Jf e 5 and X £ V* . Here V* C S* is the 'habitat' already mentioned in the 
previous subsection, a vector subspace of S* . The choice of V* should be physically motivated. 
Natural choices are V* = TCdiS or the space of vertex smooth states 29], but there may be 
others. An evident requirement here is, of course, that the choice of V* C S* should be such 
that the limit 1)5. 9|) exists. Although it is not entirely clear (to us) what precisely the conditions 
on V* should be for this limit to exist, and to what extent they would determine H*[N]x G V* 
uniquely, it seems obvious that the limit (|5.9() will not exist if V* is 'much larger' than 7 id\fj - 
There is an alternative point of view, which does not require 'habitats': according to |2^ . l74l |. 
the choice V* = Tidis in fact is singled out, and formula (|5.9() should not be viewed as defining 
a dual Hamiltonian on V* , but rather as defining the limiting Hamiltonian as a weak * limit on 
the kinematical Hilbert space Wkin w.r.t. the weak* topology defined by V* = Tidis In that 
interpretation, only the r.h.s. of (|5.9j) is meaningful — but let us repeat that there is no state 
in Wkin that could be interpreted as "lim^^o -f^[A^) e]*" -^^ 

Independently of whether or not one adopts the 'habitat' point of view, an essential ingredient 
in existing constructions is the 'diffeomorphism covariance' of the Hamiltonian. This means that 

^^We should caution readers that it is a matter of debate as to what constitutes an 'ambiguity', or whether a 
particular choice narrowing down the number of possibilities should be viewed as 'natural' or even mandatory in 
order to obtain sensible expressions (one example being the 'translation invariant' choice of measure factor in the 
definition of the volume operator). At any rate, LQG practitioners may well disagree with our terminology here. 

^^There are actually infinitely many such operator topologies w.r.t. the dual pair (<S*,<S), depending on the 
choice of V* (see e.g. 136], p. 88). The bigger the subspace V*, the stronger the topology. 

^^A more recent proposal to circumvent these difficulties is the so-called Master Constraint Programme of |2^. 
where the Hamiltonian constraint is replaced by another constraint ('master constraint') oc J d'^x (iJ)~^/^ (_ffo(x))'^ 
obtained by squaring the original Hamiltonian. This operator can be implemented as a positive self-adjoint operator 
directly on Hdiff (after suitable regularisations) , so questions about 'habitats' are sidestepped. 
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by transferring the action of H to the space of diffeomorphism invariant states, the existence of the 
above hmit can be ensured by the fact that the habitat V* is sensitive only to the diffeomorphism 
invariant characteristics of the regulator. This is also borne out by the construction of [2^ which 
takes V* to be the (larger) space of 'vertex smooth states'. (But note that from the point of view 
of the Hamiltonian action, going to the slightly larger habitat of 'vertex-smooth' states does not 
make much of a difference; namely, when acting on a spin network state, the Hamiltonian does 
not shift nodes around. Hence its action on Ti^m or on the larger habitat of 'vertex-smooth' 
states is essentially the same.) In addition to restricting the choice of V*, requiring the limit (|5.9|1 
to exist eliminates at least some of the ambiguities inherent in the definition of the Hamiltonian. 

To sum up, the Hamiltonian H[N\ can only be defined, via ()5.9() . as a weak limit. This is 
reminiscent of the LSZ formalism of ordinary quantum field theory, where one also defines the 
asymptotic field operators as weak rather than strong limits. It is remarkable that the known 
constructions make essential use of the diffeomorphism covariance, and thus implicit use of the 
constraints, already in the very definition of the Hamiltonian. However, for all we can tell, the 
question as to which V* is 'best' is not settled. On a more practical note, and independently 
of the choice of V*, we note that the concrete evaluation of the limit requires simultaneous 
knowledge of the action of the original regulated Hamiltonian H\N,€\ on all elements of S. 
Furthermore, explicit expressions for the dual Hamiltonian H*[N] that could be used for any 
practical application of the formalism are even more difficult to obtain than for the Hamiltonian 
acting on spin network states, and certainly beyond the scope of our modest efforts in this 
section. Nevertheless, we should stress once more that the dual Hamiltonian H* is our true 
object of interest, and that all further considerations (constraint algebra, semi-classical states, 
and so on) should be based on it, rather than on the regulated Hamiltonian acting on S. With 
these caveats in mind, let us now proceed. 

To motivate the form of the quantum Hamiltonian one starts with the classical expression, 
written in loop variables. Despite the simplifications brought about by equation ()3.21|) . the 
Hamiltonian constraint still looks formidable, 

H[N] = / d3xiV^^(e'^^=F^„,--(l + 72)K[„-i^„]M . (5.10) 

In order to write the constraint in terms of only holonomies and fluxes, one has to eliminate the 
inverse square root as well as the extrinsic curvature factors. This can be done [3] using the 
relations 1)3. 26() - ()3. 28(1 . Inserting these into the Hamiltonian constraint one obtains the expression 

H[N]= Jj^'xNe'^''PTr(^Frnn{Ap,V}-^il+j''){Arn,K}{An,K}{Ap,V}) . (5.11) 

This expression is the starting point for the construction of the quantum constraint operator. 

The next step is to express the Ashtekar connection Afn as well as field strength Fmn and 
the extrinsic curvature K in terms of holonomies. As in lattice gauge theory, this requires the 
introduction of a cell structure. In the classical expression one splits the target space into cells 
of size ~ in coordinate space. In order to express connections and curvatures in terms of 
holonomies, one now has to choose edges Cm within each cell, along which the holonomies are 
going to be evaluated. Once the edg chosen, one uses that the connection along a line 

of coordinate distance e can be approximated by 



I^I^A,V^ = -he[A]-^hM],V} + oie) , 
y^A,K^ = -he[A]-^{hM],K} + o{e) . 



(5.12) 
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For the field strength one has the standard lattice expression 



^ F[A] = l{hg'p^,)[A] - /iaP(.)[A]) + o{e^) , (5.13) 

where P{e) is a two-dimensional infinitesimal plaquette with boundary dP{e), that shrinks to 
zero as e — > 0. The extrinsic curvature is expressed using the relation (|3.28() . in which in one 
substitutes H5.13() . Finally, as for the volume operator, one replaces integrals with Riemann sums 
over the cells, which we will again label with a. 

In order to quantise the classical Hamiltonian H5.11() with the above regularisation, we next el- 
evate all classical objects to quantum operators as described in the foregoing sections, and replace 
the Poisson brackets in (jS-llf) by quantum commutators. The resulting regulated Hamiltonian 
then reduces to a sum over the vertices Va of the spin network with lapses N{va) 



H[N,e] = ^ iVK)6-^Tr ((/.aP_(.) - /^ai^^)) i^P'^^ 
a 

- ^(l + 7')E^(^")^""''Tr [h^,K] [hn,K] h^'[hp,V]) , (5.14) 

a 

where we have already assumed a specific ordering of the operators. The fact that the regulari- 
sation parameter e drops out 'miraculously' from this expression, and the integral ()5.11|) can be 
replaced by the sum 1)5. 14(1 . requires the Hamiltonian density to have the correct weight, and 
would not work with the weight two density (|3.21j) , cf . footnote El on page El If one does not 
want to make reference to any particular spin network, the Hamiltonian can also be expressed 
more abstractly continuous sum 

H[N,e] = Y,N{^)Hi^;e) (5.15) 

which, on any given spin network, reduces to a finite sum as in (|5.14|) . In the remainder, we 
will keep the regularisation parameter e > fixed, and not always indicate it explicitly in the 
formulas below. 

A first stumbling block here is the presence of the volume operator V in 1)5. 14() . Although its 
square has been expressed in terms of fluxes in section l4.4.21 the spectrum of this operator is only 
partially known, see, however, 41 , Ill5l . [ll7] for recent progress (recall that q is not necessarily 



positive). The further evaluation of 1)5.14(1 on a given spin network wave function would now 
in particular require a diagonal basis of spin network states, on which we could determine the 
square root of \q\ for each node, cf. (|4.3U|) . To treat all the relev ant terms, this diagonalisation 



would have to be repeated several times, along the lines of 117l |. with intermittent application 
of other operators. Although all the steps required in this calculation are thus clearly laid out, 
presently available technology does not, as far as we are aware, allow a complete evaluation of 
H[N, e] on a given, but arbitrary spin network state. For this reason, we will only be able to 
exploit certain qualitative properties of the volume operator in our analysis of the Hamiltonian 
constraint below. We will thus adopt a pragmatical approach in the remainder, performing all 
computations with the local trilinear operator q instead of the volume operator itself, and assume 
that we can divide by the relevant power of the volume at the end (thereby postponing the most 
difficult part of the calculation) Accordingly, the arguments in this section should be viewed 



on r - , „ , ^ , - , -1/2 , 



We note that replacements such as [. . . , V^] — > sgn(g) |g| [■••5(7] introduce further operator ordering 
ambiguities into the definition of the Hamiltonian operator. 
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as no more than a gedanken calculation, and there is no claim that this is the way the calculation 
should ultimately be done. 

Apart from obvious ordering ambiguities, and the difficulties with volume operator just de- 
scribed, the action of the quantum expression 1)5.14(1 on a generic spin network state also depends 
strongly on the choice of the coordinate system on S. To be more precise, the action depends 
on whether or not the basis vectors {i.e. the directions m,n and p) are aligned with the spin 
network edges at a given node. There exists a large number of possibilities for these alignments, 
which are not diffeomorphically equivalent, and one might thus simply define the action of the 
Hamiltonian by employing a particular alignment, as LQG does. Some of these possibilities are 
explicitly listed below. We regard these choices as genuine ambiguities in the quantisation pro- 
cedure, in the sense that observables quantities may depend on them, if they can be shown to 
exist at all. This potential prescription dependence, and the abundance of possible outcomes for 
physical quantities are most difficult to digest from the perspective of conventional quantum field 
theory, especially because, the result for the corresponding classical quantity does not depend on 
such choices when the limit e — > is taken. 

To get an idea of the complications, let us consider the action of the first term in the regulated 
Hamiltonian constraint, which is occasionally called the 'Euclidean Hamiltonian' (because with 
a Euclidean signature and 7 = ±1, the bothersome second term in (|5.14|) would be absent, as it is 
for 7 = ibi with Lorentzian signature). We will act with this operator on a three- valent vertex v 
in a spin network, for given e > 0. The spin network wave function on which we will act is given 

by 

^ = m^.,, ma.,, ma,,, • • • - (^-le) 

where the dots denote other parts of the spin network. Since the Hamiltonian acts node by node, 
we can restrict our attention to this local structure at v. Let us also write out all indices on the 
Hamiltonian operator, restricted to that part which acts at v, 

H = sgniq)\r'/'e-^FU^^Jhr)-' \q, {hi) + . . . , (5.17) 



where the dots now denote terms which act at other vertices. The choice of the spin j is a 
new ambiguity. In the first step, the operator h~^[hp, q] = q — h"^ qhp acts at the vertex of the 
network. When q in the second term acts on three holonomies already present in the spin network, 
the result cancels against the first term. (Note that the result only vanishes because we are using q 
instead of V, so that the local volume operator acts tri- linearly in holonomy derivatives. This is 
no longer true when one uses the square-root expression for V in the Hamiltonian constraint). 

When q acts on two existing holonomies and on the hp factor, however, the result is non- 
trivial. In order to work out this non-trivial action of q, it is necessary to fix the orientation of 
the coordinate system on S. We will here choose the m, n and p directions to be aligned with 
the three edges 61,62 and 63 of the spin network respectively, as is usually done in the LQG 
literature. Note, however, that this is not the only possibility; we will return to this arbitrariness 
shortly. In effect, this means that the quantum Hamiltonian is postulated to be 



a I,J,K 



^ = E E ^(^") ^""^ Tr [{hep,, - hg'pj h;^ [hs,,V] 

J,K 

^(1 +7')E E N{va)L''''TT (/i^; [K„K] h-} [hs,,K] h-l[K,,,V]) . (5.18) 



a I,J,K 



^ AxiA thus often designated as 
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The second sum is a sum over all triplets of edges emanating from the vertex Va in the centre of 
the cell a. The boundary of a plaquette in the plane spanned by the edges sj and sj is denoted 
by dPjj. Finally, the value of the constant, fully anti-symmetric tensor L^^^ depends on the 
particular vertex a. 

Having made this coordinate choice, we can now return to our calculation. The situation is 
as in the first graph of figure IIUI Because the small edge p is extending in the 63 direction, we 
can use the tensor product decomposition rule to first 'merge' the hp holonomy with (a segment 
of) the /ie3 holonomy of the spin network. For two arbitrary holonomies along identical edges, 
this decomposition rule reads 

(^e')ai/3i (^e')a2/32 = X] ^"i "2^1 '^£/3272 ^^^^) ^1-12 ' ^^'^^^ 
k 

By splitting the edge 63 into a 'small' piece, overlapping with the p edge, and a 'larger' 
remainder, we can write 

where the edge 63 is the small piece of the edge, equal to p. Applying the tensor product rule 
to the holonomies hp and /ig^ one now obtains (recalling that also m and n are aligned with ei 
and 62 respectively) 



Hi. = sgn(g) i^rV2 {FU,,,M^;l,{h^^,,j 



'3' "373 

X]^A373P3^Aifto-3 (^e|)p3a3 (^eOaift (^62)02/32 '^^iftjft ' ^^■'^^'> 

The factor in brackets corresponds to the two merged holonomies. We can now work out the 
action of the first line of this expression. 

As explained in section 14.4.21 the square of the local volume operator inserts a r matrix at 
the end of the three holonomies on which it acts. In our case, this means that the matrices 
are connected to the indices Ai, /3i and /32, as these indices correspond to the end-points of the 
three holonomies which are located at the vertex v (remember that in expression (t5T7|) . the 
only non-zero contribution arises when q acts on two holonomies of the spin network and on 
the holonomy (hp)\.^\^ inside the Hamiltonian constraint; the Ai index is the one located at the 
vertex). Up to a normalisation factor, this leads to 



iii, = sgn(g) e"^^ (i^^^J^^^^ {h{,^Y^^^^ {h^)^^^^ 



^''). \YCi^l^^&^^„^UhH) (hi]T^) AKl^") « (^^Tfl •••• (5-22) 

/A1T3 1 A373P3 r3ft(T3 V e'.U pzcTz I V /«i/3iV £2 /a2/32 P1P2P3 ^ ' 



63 ^P30"3 



What is left to do is to insert the plaquette. Because of our choice of coordinate system on S, 
the plaquette extends in the plane spanned by the ei and 62 edges. In the limit e ^ 0, the field 
strength can be expressed in terms of elementary holonomies as 

(-^eie2)AiA2 ~ 



We put quotation marks here, because there is no a priori reference metric w.r.t. which the lengths of these 
subedges could be measured. 
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where /ieie2 denotes a holonomy in the direction which connects the edge ei with the edge 62- 
When we insert this expression in ()5.22|) . we can again spht the edges in a small and large piece 
as in (|5.2U|) and then use the tensor product decomposition H5.19|) to merge the two holonomies 
in the ei direction and the two holonomies in the 62 direction. This introduces four new Clebsch- 
Gordan coefficients, with the result 

(^4')a2A3 (^eie2)KiK2 '^nrafti (^'')/3iri (^^) /32r2 (^") Aira X] '^AgTsPS (^e|)p3fT3 



X 



/^j jlkl i^j jlki (uki\ 1 I /^j j2k2 /^j j2k2 (uk2\ 1 (^'}A\ 

Z^'-^Ai/3io-i'^Kl7lpA"'eiVpi<7i I '-'X2P2<T2^ >^2l2P2\'^ep p^aij ■■■■ ^^■^^> 

Altogether, these manipulations have introduced three new vertices, at locations e away along 
each edge, and they have also modified the original vertex. The three new vertices are given by 

('T/(i)'\fciiii = rjjjiki 

/'T/(2h'^2i2j _ nj j2k2 (r oc\ 

! 92121^2 ~ ^K272P2 ' yO.^Oj 

(y{'i)\k3.m _ (-ij j'iks 

while the old vertex has been changed to 

^fclA:2fc3 i ^(^ijlfcl (^ij2fc2 ^ii3fc3 ^ a6c/6N /cN /aN C'^j^^,' . (5.26) 

cri(T2CT3A2 Al/3i(Ti A2/520-2 Ta/Jsas \ / [3iti \ //32T2^ ''A1T3 T1T2/33 V ' 

One can still merge {h-^^") x^x^ ^^'^ ^^ei') p^a^ ^^i^S tensor product decomposition rule, which 
will reduce the vertex V^^'^ to an overall factor times a Kronecker delta. We will refrain from 
spelling out these details. 

Let us now comment on other possible coordinate systems and alignments. The use of the 
tensor product decomposition 1)5. 19() depends crucially on the edges of the spin network being 
aligned with the edges of the plaquette (and, through (|5.14|) . with the coordinate system on S). 
If this is not the case for the m,n and ei,e2 edges, the plaquette will not get glued into the spin 
network, but instead remains freely floating, only connected at the original vertex (see figure [TT|) . 
Similarly, if the p edge is not aligned with 63, the plaquette will only be connected to the spin 
network through the matrix in the construction above. In the LQG literature, such cases are 
excluded by hand, by always choosing a local coordinate system around each vertex which is 
completely aligned with (three of the) edges emanating from the vertex. This choice is usually 
justified by invoking background independence, which would be violated otherwise 

We have so far not discussed the second term in the constraint (|5.14l) . which arose from the 
product of extrinsic curvatures in (|5.1fl|l . This term is hardly ever discussed in detail, and with 
good reason: it is far more complicated than the first one. In fact, the technical complications 
which we have just described pile up in such a way as to make any concrete calculation almost 
unfeasible. The complexity of this term stems from the complexity of the intermediate expression 

[Ar,K] = h-} [/ie,,e^^^ [{K.,^ - h-]jh~^ [he„V], v]] , (5.27) 

which is obtained using ()3.26|) - H3.28() . One recognises the combination of operators which opens 
up the vertex and inserts a plaquette, but there there is an additional volume operator which 
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Figure 10: The action of the first term in the Hamiltonian constraint on a spin network, in 
the special case that the m, n and p directions are ahgned with three existing edges of the spin 
network. The three edges are drawn in the plane for convenience, but should span a tetrahedron 
in order for the volume operator to act non-trivially. Filled dots denote insertions of r matrices. 
In the last two steps, we have only drawn one term in the plaquette H5.23|l . 

applies further changes to the node. If the coordinate system on S is aligned in such a way that 
all three directions coincide with an existing edge, the action of the second term in the constraint 
is roughly to glue in three new edges, introducing six new vertices. But as the discussion above 
shows, there is a plethora of other possibilities. 

Let us, for clarity, summarise the problems which we see with the Hamiltonian constraint. 

• Ambiguities. The construction of the Hamiltonian outlined here is plagued with a variety 
of ambiguities. Firstly, a choice has to be made for the coordinate system on S. The 
action of the Hamiltonian depends in a qualitative way on this choice. Secondly, there 
is an ambiguity in choosing the representation in which the trace is evaluated [i.e. the 
representation in which the field strength Fmn and the holonomy hp are taken). Thirdly, 
there are operator ordering ambiguities. Finally, and less conspicuously, there may exist h- 
dependent 'higher order' corrections to the Hamiltonian (|5.14|) . which are still compatible 
with diffeomorphism covariance, but would not affect the classical limit. Not only will 
physical predictions depend on the choices made to fix these ambiguities, but these choices 
are important in order to determine whether or not the limit e ^ of 1)5. 9(1 exists at all. 

• Ad-hoc choices. The quantum constraint should be well-defined and unambiguous on the 
full space of spin-network states. If this is not the case, it should be shown that the choices 
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Figure 11: The action of the first term in the Hamiltonian constraint, for the case in which 
the m and n directions on S are not aUgned with any of the edges emanating from the vertex. 



do not matter for the physical predictions, or alternatively that a particular choice is singled 
out by physical consistency. By focussing only on the operator (|5.18j) , it becomes impossible 
to tell whether the constraint operator is indeed independent of the coordinate choice on S. 

Ultralocality. From the explicit discussion of the Hamiltonian constraint it should be clear 
that its action is always 'ultralocaP, in the sense that all changes to the spin network are 
made in an e — > n eighbourhood of a given vertex, while the spin network graph is kept 
fixed 137 . 138l . 139l |. More specifically, it has been argued that the Hamiltonian acts 



at a particular vertex only by changing the intertwiners at that vertex, viz. 

H{v) |r, {i}, {C}) = f^lr, {j}, Ca{T, {j}, I, J)) (5.28) 

where [/" is the traceless part of the holonomy associated with the plaquette P/j and Ca 
a triplet of new intertwiners at u , and the sum is over pairs of edges /, J entering the 
vertex. This is in contrast to what happens in lattice field theories. There the action 
of the Hamiltonian always non-trivially links two existing nodes, the plaquettes are by 
construction always spanned between existing nodes, and the continuum limit involves the 
lattice as a whole, not only certain sub-plaquettes. 



6 Constraint algebra 

As we already emphasised in the introduction, the proper closure of the quantum constraint 
algebra is perhaps the crucial consistency requirement that must be met by any candidate theory 
of quantum gravity. Of course, the devil hides in the meaning of the word 'proper'. To find out 
where it might possibly hide, we now return to the points raised at the end of section We will 
not bother here with the commutators involving the Gauss constraints, which are straightforward 
to take care of. The canonical (Poisson) brackets involving the diffeomorphism generators read 

{Dra[M] , D,,[N]\ = I d^X {M^dmN'' - N'^d^M'')Dn , (6.1) 
{d^[M] , Ho[N]j^^ = jj?x {M'^dmN - NdmM^)HQ . (6.2) 
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Here = ^m°'Da is the diffeomorphism generator, while the smeared Hamiltonian and smeared 
difFeomorphism generators are defined as 

Ho[M] := [ d^xHoix)M{x) , D^[M] := [ d^x L>„(x)M™(x) . (6.3) 

Because infinitesimal difi'eomorphisms do not actually exist on Wkinj only the exponentiated 
version of (|6.1|) makes sense in LQG. For the same reason, relation (|6.2j) is implemented by 
considering only finite (exponentiated) diffeomorphisms (p. Given a realisation of 99 on the Hilbert 
space through a unitary operator U[ip], the bracket (|6.2|) can be exponentiated formally to^^ 

U[^] Ho[N] U[ip]-^ = Ho[ip*N] . (6.5) 

These relations constitute what we would call the 'kinematical part' of the constraint algebra. 

The crucial, remaining relation is the commutator of two Hamiltonian constraints, whose 
classical analog reads 

\Ho[M] , Ho[N]} = [ d^x {Md^N - Nd^M)g"'''Dn • (6.6) 

Together with the 'kinematical relations' above, this relation encapsulates the closure of the 
classical algebra of spacetime diffeomorphisms in the Hamiltonian formulation (we should note, 
however, that the relation between timelike diffeomorphisms on the one hand, and the 'motion' 
generated by the Hamiltonian constraint on the other, is quite subtle J141]). To be sure, the 
proper closure of this part of the quantum constraint algebra, i.e. the issue of promoting (|6.6jl 
to a relation between operators acting on a suitable Hilbert space of states, is also one of the 
major unsolved problems of the conventional geometro dynamics approach. There it has not been 
possible to resolve the ordering and other ambiguities, and to get rid of the infinities resulting 
from the clash of functional differential operators (the situation is partly better, though, in certain 
supergravity models whose constraint superalgebra always contains the bosonic constraint algebra 
as a subalgebra ^76]). 

6.1 The bosonic string as an example 

Before continuing our discussion of the quantum constraint algebra in canonical gravity, we would 
like to recall that there does exist a model which can be solved exactly, and all the way. Hence 
it can serve as an example of how quantisation should work. This model is the bosonic string 
(in a Polyakov type formulation), wh en v iewed as a bona fide model of matter coupled quantum 
gravity in two dimensions (see also |l42l | for an earlier discussion of the bosonic string in this 
context). For this purpose, we only need to forget about the target space interpretation and its 
ramifications, such as string scattering, or higher order corrections in involving splitting and 
joining of strings (in the present context, the latter processes would amount to changes of the 
spatial topology, 'baby universes' and the like, hence, 'third quantisation'). The resulting theory 
is manifestly independent of the two-dimensional space-time background, free of divergences and 



■^'^When this relation is understood to act on a state in H diff, one can use invariance of the state under diffeo- 
morphisms, tf[tfi\\'4>) = \il>), to rewrite equation H6.5|l as |l4fi | 



[U[^] , 7^0 [iV]] = Mf'N " N\ . (6.4) 
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other pathologies, and also possesses nice observables (vertex operato rs). Mor e sp ecifically, the 
diffeomorphism and scalar constraints D and H are given by (see e.g. |l43| | or [l4J]) 



D[f] = \j\af{a)P^{a)d,X^^{a), 

H[N\ = ^l\aN{a)(^P^{a)P>'{a) + d^X^{a)d^X''{a)) . 



(6.7) 



with shift /(cr) and lapse N{a). The classical constraint algebra, given by (|6.1|) . (|6.2|) and (|6.6|1 
in four dimensions, reduces to 

{D[h], D[h]]^^ = \D[hdh - f2dh] , 

{D[f], H[N]}^^ = \H[fdN - Ndf] , (6.8) 
{h[Ni], H[N2]]^_^ = ^D[NidN2 - N2dNi] . 

Canonical quantisation is straightforward in the conformal gauge. The state space of the theory 
is simply the tensor product of the space of one-particle wave functions and an infinite number of 
harmonic oscillator Hilbert spaces associated with the excited string modes. Composite operators 
become well defined, i.e. have finite matrix elements between arbitrary states in J^, through 
normal ordering. Both the kinematical and the Hamiltonian constraint are hermitean operators. 
With the normal ordered energy momentum tensors T±±(a) = : {H{a) it D{a)) :, the classical 
constraint algebra is modified by the anomaly (the central term), which is responsible for almost 
everything that is non-trivial about string theory, 

[T±±(a) , r±±(a')] = S'{a,a'){Tii{a) + r±±(a')) + hc6"'{a,a') . (6.9) 

Our main point in repeating this well known story here is that this algebra holds on the full 
space of states, i.e. prior to the imposition of any constraints, hence closes off-shell. Only at 
this point can we impose and solve the quantum constraints to obtain the subspace J-q C oi 
physical states — with the well known result ('No Ghost Theorem') that negative norm states 
are absent if the number of scalar fields is < 26. Because the Hilbert spaces J-'q and J-" of physical 
and unphysical states, respectively, are well defined, there is no need to discuss 'habitats' here. 
Of course, one key simplifying feature of the string model must be pointed out here, which is 
not shared by (matter coupled) gravity in higher dimensions: (|6.8j) is an honest (albeit infinite 
dimensional) Lie algebra (and ()6.9|) closes in Dirac's sense), whereas (|6.1() . 1)6. 2|) . and 1)6. 6() are 
not. This Lie algebra furthermore possesses a triangular decomposition, allowing for a systematic 
study of its representations, and this fact is a crucial ingredient in establishing the above results. 

Very recently, LQG methods have been applied to string theory, but have been s hown to l ead 
to results which are at variance with those obtained using Fock-space quantisation 145l . 146l |. It 



therefore does not appear that LQG can recover conventional results; instead it gives rise to an 
inequivalent quantisation of the string model, in which the quantum constraint algebra is realised 
without a central term, and without any restriction on the ta rget space dimension. The meaning 



of these results is currently subject of an intense debate, see |147l . Il48l |. 

Another simple model which can be quantised all the way is dilaton gravity in 1-1-1 dimensions. 
For this model the constraint algebra has been worked out for various inequivalent representa tions 
of the op erators, an d the physical consequences have been analysed in detail. See for instance |l4g| . 
15ril . 151 1 . See also |l52l | for a review of dilaton gravity with many further references. 
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6.2 On-shell vs. off-shell closure 

We have seen that LQG treats the diffeomorphism and Hamiltonian constraints in a very different 
manner. Diffeomorphism invariance is implemented by an averaging procedure that makes the 
states invariant under finite (i.e. 'exponentiated') diffeomorphisms. In LQG a canonical gener- 
ator of diffeomorphisms simply does not exist, since the lack of weak continuity that goes with 
the scalar product (|4.7|) prevents us from differentiating finite diffeomorphisms so as to obtain 
'infinitesimal' ones. By contrast, at least so far, group averaging is of no use, even at the formal 
level and allowing for whatever approximation scheme, when it comes to imposing the Hamilto- 
nian constraint. Therefore the latter, whose cumbersome form we have exhibited in the foregoing 
section, must be analysed 'the traditional way' This restriction applies also to the commuta- 
tor of two Hamiltonian constraint operators. A further (and well known) difficulty is that, even 
for the classical theory, this commutator does not generate an ordinary Lie algebra, in the sense 
that there appear field dependent structure 'constants'; as is well known, these lead to ordering 
problems and ambiguities in the quantum constraint algebra. 

While there is general agreement as to what one means when one speaks of 'closure of the 
constraint algebra' in classical gravity (or any other classical constrained system), this notion is 
more subtle in the quantised theory. Let us therefore clarify first the various notions of closure 
that can arise: we see at least three different possibilities. As we explained in the foregoing 
section, the 'true' Hamiltonian H* is defined as a weak limit on some habitat V* C S* by means 
of (|5.9j) . Assuming that this subspace is preserved by the action of H* , the strongest notion is 
'off-shell closure' (or 'strong closure'), where one seeks to calculate 

[H*m , H*[N2]] = d*{Nr,N2) (6.10) 

without further restrictions on the states, on which this relation is supposed to hold — in other 
words on a sufficiently large 'habitat' V* that does not only contain diffeomorphism invariant 
states. The operator 0*{Ni;N2) here is the quantum analogue of the r.h.s. of (|6.6|) . and would 
be proportional to the diffeomorphism generator (but this commutator might also differ from 
the classical one by certain quantum modifications). Together with the 'kinematical part' of the 
quantum constraint algebra, the above relation would express the quantum space-time covariance 
of the theory. Strong closure is realised for the string model, with the well known result ()6.9() . 
If one could resolve the ordering ambiguities and eliminate all singularities, this is also the type 
of closure expected to work in conventional geometrodynamics. On the other hand, although 
the calculation in principle does make sense in LQG, too, the r.h.s. of (|6.1fl|) may not exist as 
an operator in LQG, as we explained above. Of course, even if the diffeomorphism generator 
does not make sense as an operator, it is still conceivable that a combination like g^^^Hn does 
exist as a 'composite operator'. In this case, however, this operator would not generate a closed 
algebra: repeated commutation would generate an infinite tower of new operators, and one would 
therefore no longer be discussing the algebra of constraints. 

At any rate, it appears that the goal of determining 0*(A'^i; A'^2) as a bona fide 'off-shell' 
operator on V*, and prior to the imposition of any constraints, is unattainable within the current 
framework of LQG. For this reason, LQG must resort to weaker notions of closure, by making 
partial use of the constraints. More specifically, equation can be relaxed substantially by 

demanding only 

[H*[Ni], H*[N2]]X = 0. (6.11) 



To overcome the difficulties, it has even been suggested that in analogy with the diffeomorphism generator, 
the Hamiltonian itself may actually not exist as an infinitesimal generator, but that only 'finite' translations in 
time might be well defined [S^, il53| . 
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This 'weak closure' should hold for all states X in a restricted habitat V* of states that are 
'naturally' expected to be annihilated by the r.h.s. of 1)6. 1U|) . and that are subject to the further 
requirement that the Hamiltonian can be applied twice without leaving the 'habitat'. The latter 
condition is, for instance, met by the 'vertex smooth' states of "2^. As shown in |2^. Ei!]. the 
commutator of two Hamiltonians indeed vanishes on this 'habitat', and one is therefore led to 
conclude that the full constraint algebra closes 'without anomalies'. 

The same conclusion was already arrived at in an earlier computation of the constraint algebra 
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23, [TJ], which was done from a different perspective (no 'habitats'), and is based on the choice 
V* = Ti-diS, the 'natural' kernel of the r.h.s. of (|6.10|) . Here the first step consists in calculating 
the commutator of two regulated Hamiltonians (see section 

d{Nv,N2; e) := [H[Ni,e] , H[N2, e]] , (6.12) 

for fixed but arbitrary e > on the space S of spin network states. As for the Hamiltonian itself, 
letting e — > in this commutator produces an uncountable sequence 0(A'^i, A''2; e)^' of mutually 
orthogonal states w.r.t. the scalar product (|4.7j) : consequently, lime-^o O(-^i) -^2; e) again does 
not exist in the usual sense, but only as a weak* limit. More specifically, one shows that 27, 3| 



lim (A:'|0(iVi;A^2;e)^) = (6.13) 

for all X G Wdifr, and for all ^' € 5 (the 'diffeomorphism covariance' of the Hamiltonian is 
again essential for this result); in this sense, the constraint algebra is free of anomalies. Let us 
stress that ()6.11() and ()6.13|) are by no means the same: in 1)6. 11() one uses the dual Hamiltonian 
defined in (|5.9j) (where the limit e ^ has already been taken) and performs the calculation on a 
subspace V* C 5*, whereas the calculation of the commutator in (|6.13() takes place on the space 
S, that is, inside TCkim and the limit e ^ is taken only after computing the commutator of 
two regulated Hamiltonians on S. In other words, these two operations (taking the limit e — > 0, 
and calculating the commutator) need not commute. Because with both 1)6.11(1 and (|6.13() . one 
forgoes the aim of finding an operatorial expression for the commutator [ff*[A'^i], i?*[A''2]] as an 
operator on V*, or the commutator (|6.12() on the kinematical Hilbert space, making partial use 
of the constraints, we say (in a partly supergravity inspired terminology) that the algebra closes 
'on-shell'. 

At this point we should like to remark that computations done so far make only little use 
of the detailed structure of the terms in the Hamiltonian. In particular, (on-shell) closure is 
achieved independently of whether or not one takes into account the second term in (|5.14j) . 
and there may exist (diffeomorphism covariant) modifications of the type discussed at the end 
of this section, or 'higher order' modifications of the Hamiltonian ()5.14|) . which do not alter 
this result. The analysis of j^l uses only schematic formulas like (|5.28|) rather than explicit 
operatorial expressions for the Hamiltonian. All this is contrary to naive expectations, because 
computations of more complicated gravity-matter systems (such as the supergravity constraint 
algebra studied in [^) show that the detailed structure and the interplay of all terms in the 
constraints is indispensable already for the proper closure of the classical constraint algebra. 
This may indicate that crucial information is lost, and that, in a sense, the on-shell closure is 
already 'built into' the formalism via the diffeomorphism covariance of the Hamiltonian, and by 
restricting the calculation to special states (or matrix elements) only. Still, one would hope to 
be able in principle to extract from this computation some kind of explicit operatorial expression 
for the commutator. However, (at least our) attempts to actually evaluate ()6.12|) on a given 
spin network wave function ^ along the lines of section 15.21 and before transferring the result 
to V*, get stuck quickly: the resulting state would be substantially more complicated than the 
(incomplete) expressions displayed there. 
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The procedure of making partial use of the constraints prior to the evaluation of the constraint 
algebra, as well as in the very construction of the Hamiltonian constraint, is very different from 
the conventional treatment of quantised canonical systems fSQi], as exemplified by the bosonic 
string in the previous subsection. There, one first studies the constraint operators and tries to 
make them well defined as quantum operators without constraining the allowed states, and before 
checking the closure of the constraint algebra. When one then in a second step computes the 
algebra, one usually encounters subtle quantum (i.e. ^-dependent) modifications of this algebra, 
whose determination requires great care. Because these modifications are subject to certain 
constraints (Wess-Zumino type consistency conditions and their offspring in BRST cohomology), 
they can often be determined in a representation and gauge independent way, as is the case for 
the group Diff(S'^). If one succeeds in working out this algebra and demonstrating its consistency, 
the quantum constraint algebra closes 'off-shell'. Only after ensuring the closure of the quantum 
constraint algebra does one proceed to impose and solve the quantum constraints. Let us also 
remind readers that, at least according to conventional (quantum field theory) wisdom, 'on-shell' 
calculations of the algebra might be 'empty': in general it is not co rrect to use a symmetry before 
it has been shown that it can be implemented without anomalies |l54l |. 

Why do we emphasise off-shell closure, despite the difficulties (ordering ambiguities, field de- 
pendent structure constants, and the like), which likewise plague other approaches to canonical 
quantum gravity? Our main reason for emphasising this requirement is two-fold: first and most 
importantly, the off-shell closure of the quantum constraint algebra (which, as we pointed out, 
need not coincide with the classical algebra) encapsulates the key property of quantum gravity, 
namely quantum space-time covariance. Secondly, this property furnishes an excellent means to 
distinguish between the 'correct' theory and a mere regularisation: consider for instance replacing 
the kinetic part of the WDW operator by a smeared point-split operator. With this modification, 
all constraint operators become formally well defined, so one might contemplate 'defining' the 
quantised theory in terms of them, and simply proceed to study the solutions of these modified 
constraints. This procedure would, of course, suffer from the same kind of ambiguities encoun- 
tered in section [2121 and before, but it might still take us a long way if we did not pay due attention 
to the constraint algebra. If, however, we do subject the modified WDW Hamiltonian to this 
acid test, it is immediately clear that the resulting operator (or any other obtained by such a 
regularisation) is 'wrong': it does not satisfy the correct algebra, whence the quantised theory 
would not satisfy the basic requirement of quantum space-time covariance. 

As another example, consider modifying the Hamiltonian constraint of string theory by mul- 
tiplying it with an operator which commutes with all Virasoro generators, 

(f++ + f__) {f+++f^^)C, (6.14) 

where [C ^Tj^j^] = and [C',T ] = 0, while keeping the diffeomorphism constraint (2++ — T^^) 

unchanged. There are many such operators C in string theory. The simplest is to take C to be 
equal to the mass operator P^P^ minus some arbitrary positive integer, where P^ is the conserved 
global momentum of the string (which indeed commutes with all Virasoro operators). In this 
way, we arrive at a realisation of the constraint operators which is very similar to the one used in 
LQG: the algebra of spatial diffeomorphisms is realised via a (projective) unitary representation, 
and the Hamiltonian constraint transforms covariantly as in (|6.5|) (the extra factor C does not 
matter, because it commutes with all constraints). In a first step, one can restrict attention to 
the subspace of states annihilated by the diffeomorphism constraint, the analog of the space TLdm- 
Imposing (T++ + T ) Clphys) = on this subspace would now produce a 'non-standard' spec- 
trum by allowing extra diffeomorphism invariant states of a certain prescribed mass, but the 
algebra would still close on-shell, i.e. on the 'habitat' of states annihilated by the diffeomorphism 
constraint. The point here is not so much whether this new spectrum is 'right' or 'wrong', but 
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rather that in allowing such modifications which are compatible with on-sheh closure of the con- 
straint algebra, we introduce an infinite ambiguity and arbitrariness into the definition of the 
physical states. In other words, if we only demand on-shell closure as in LQG, there is no way 
of telling whether or not the vanishing of a commutator is merely accidental, that is, not really 
due to the diffeomorphism invariance of the state, but caused by some other circumstance. 

This, then, is our main point: by weakening the requirements on the constraint algebra and 
by no longer insisting on off-shell closure, crucial information gets lost. This loss of information 
is reflected in the ambiguities inherent in the construction of the LQG Hamiltonian. It is quite 
possible that the LQG Hamiltonian admits many further modifications on top of the ones we have 
already discussed, for which the commutator continues to vanish on a suitably restricted habitat 
of states — in which case neither (|6.11|1 nor (|6.13|1 would amount to much of a consistency test. 

7 Conclusions 

String theory and LQG pursue the same goal, a consistent theory of quantum gravity, though with 
very different means. Both approaches address core issues of quantum gravity, but concentrate 
on complementary aspects of the problem, and have led to valuable insights. For this reason, 
the opinion has been voiced that string theory and LQG may ulti matelv m erge together in a 
grand synthesis, or that LQG might become part of string theory |l55l . Il5fi| |. On the basis 
of the available evidence, this does not appear a likely outcome to us: the basic hypotheses 
underlying string theory and LQG seem impossible to reconcile without major modifications in 
either approach. There is a basic clash here, as reflected in such questions as to whether or not 
supersymmetry is necessary for a consistent quantisation of gravity. One might thus say that the 
differences between the two approaches are such that, in the end, at most one 'can be right'. 

We have reported here on the status of loop quantum gravity from a somewhat uncommon 
perspective. Our general conclusion is that, despite the optimism prevalent in many other reviews, 
more attention should be paid to basic aspects and unresolved problems of the theory. As many 
LQG practitioners may find this assessment too harsh, let us therefore summarise once more the 
basic point we have been trying to make in this review. The main issue here is not so much 
specific details of the LQG approach, where there has been a lot of serious work and considerable 
progress, but rather the question: what does one mean when one speaks of a consistent theory of 
quantum gravity? And what are the basic properties that such a theory should satisfy? It is here 
that opinions differ. In the introduction, we have spelled out some of the criteria that we believe 
such a theory should meet. Apart from the question of the semi-classical limit (which is generally 
acknowledged to be a main outstanding problem of the LQG approach), and our insistence on 
the key role of the perturbative two- loop divergence (|l.lj) . we have emphasised the importance of 
finding a good criterion for ensuring the space-time covariance of the theory at the quantum level 
— after all, this, rather than merely covariance under spatial diffeomorphisms, is the essence 
of Einstein's theory. We have recalled the well known fact that, in the canonical approach, 
it is the constraint algebra that encodes the information about space-time covariance, and that 
everything hinges on how this algebra is to be implemented in the quantum theory. Whereas LQG 
proponents generally seem to be content with ()6.11|) or 1)6. 13() as an expression of covariance, as 
well as the procedure of solving constraints 'in steps', using the diffeomorphism constraint already 
in the very definition of the Hamiltonian, we have tried to explain why we consider both (|6.11l) 



Finding a connection between classical two-form gravity in four dimensions and topological string theory |157T | 
can hardly be considered evidence for the connection between LQG and string theory. String theory provides a 
large variety of topological theories, of which the ones that appear in LQG are only special examples. More 
importantly, the quantisation methods employed there are entirely different from those of LQG. 
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and (|6.13|) as too weak, hence insufficient, in this regard. Instead, we have proposed the off-shell 
closure of the quantum constraint algebra as a criterion of space-time covariance. We have also 
argued that imposing off-shell closure may further reduce the large number of ambiguities present 
in the formalism, and in the definition of the Hamiltonian constraint, in particular. As we have 
stressed repeatedly, it is this constraint which lies at the core of canonical quantum gravity. 

In summary, we feel that there are still important problems at a basic level that need to be 
addressed and resolved before one can tell whether or not the loop quantum gravity programme 
can succeed. We hope that the present paper will provide an incentive to re-focus attention to 
these basic issues, and contribute to the debate between the different approaches. 
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